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ABSTRACT 


This  report  describes  the  design  and  evaluation  of  seismic 
classifiers  for  distinguishing  among  humans,  heavy  trucks, 
armored  personnel  carriers,  helicopters,  and  C-131  aircraft. 
The  data  used  to  develop  these  classifiers  consisted  of  many 
digitized  seismometer  responses  to  each  of  the  Intrusion  tar¬ 
gets  and  was  collected  by  the  Sensor  Development  Section  of 
the  Surveillance  and  Control  Division  (DCTI)  at  the  West  Lee 
Test  Site.  The  Interactive  Processing  Section  of  the  Informa¬ 
tion  Sciences  Division  (ISCP)  analyzed  this  waveform  data  and 
extracted  an  Initial  set  of  48  features.  The  On-Line  Pattern 
Analysis  and  Recognition  System  (OLPARS)  was  then  used  to  de¬ 
velop  several  seismic  classifier  designs  which  are  based  on 
different  subsets  of  the  initial  48  features. 
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SECTION  I 
INTRODUCTION 


This  report  documents  the  first  attempt  to  design  linear 
classification  logic  based  on  seismic  waveform  data  collected 
at  RADC's  West  Lee  Test  Site.  This  decision  logic  was  designed 
to  distinguish  among  humans,  heavy  trucks,  armored  personnel 
carriers,  helicopters,  and  C-131  aircraft.  The  classifier 
design  procedure  employed  the  following  sequence  of  tasks: 

.  Data  Collection 

.  Development  of  an  Interactive  Graphics  Tool  for 
Data  Analysis 
.  Data  Analysis 

.  Development  of  an  Automatic  Segmentation  Algorithm 
.  Feature  Hypothesis 
.  Feature  Extraction 
.  Feature  Evaluation 
.  Classification  Logic  Design 

.  Testing  Classification  Logic  with  Independent 
Test  Data 

This  effort  was  conducted  in  support  of  Project  692B 
of  the  Advanced  Sensor  Development  Program.  The  seismic 
classifier  designs  discussed  in  this  report  are  based  entirely 
on  data  collected  by  the  Sensor  Development  Section  at  its  West 
Lee  Test  Site  and  made  available  to  the  Interactive  Processing 
Section  for  completion  of  the  remaining  tasks  of  the  classi¬ 


fier  design  procedure. 


SECTION  II 

DATA  COLLECTION  AND  CONVERSION 


The  seismic  data  base  used  to  design  the  October  1st 
Instrusion  Classifier  was  collected  and  digitized  by  DCTI  at 
the  West  Lee  Test  Site  in  the  first  quarter  of  1972.  The 
five  classes  of  intrusions  observed  were  helicopters,  armored 
personnel  carriers,  C-131  aircraft,  heavy  trucks,  and  humans. 
Background  data  was  also  collected  and  used  to  determine  a 
suitable  detection  threshold  for  the  turn-on  criteria. 

The  test  procedure,  following  calibration,  consisted  of 
running  the  Intruders  along  one  of  five  specific  paths  at 
several  known  constant  speeds.  Each  intrusion  was  repeated 
with  each  speed,  path,  and  direction  as  a  check  on  repeat¬ 
ability.  Each  intrusion  involved  only  one  object,  with  the 
exception  of  humans  where  there  were  multiple  as  well  as 
single  intrusions.  Each  sensor  was  a  three-axis  low  frequency 
geophone,  Geo-Space  Model  VLF-LP-3D,  with  one  vertical  and 
two  horizontal  axes  (parallel  and  perpendicular  to  the  intru¬ 
sion  paths) . 

As  the  intrusions  were  taking  place,  the  seismic  trans¬ 
ducer  signals  were  relayed  by  underground  cables  to  the  site 
control  center,  digitized,  and  recorded  directly  on  digital 
tapes  in  the  BAMKI  format.  This  format  is  capable  of  packing 
45  simultaneous  sensor  waveforms  on  the  tape.  Since  the 
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experiment  involved  3  three-axis  geophones,  the  BAMKI  format 
sparsely  packed  9  of  1-5  possible  channels  of  digitized  wave¬ 
forms  on  each  tape.  Each  file  on  the  tape  contained  one  intru¬ 
sion  run,  consisting  of  a  number  of  320  sample  records.  First, 
the  data  was  filtered  at  500  Hz  and  then  digitized  at  1000 
samples  per  second.  Each  sample  value  was  quantized  to  any  of 
1023  values  ranging  from  -2044  to  +2044  in  steps  of  4.  The 
corresponding  strip  chart  range  was  +  10  volts  maximum.  Each 
test  usually  produced  three  or  four  magnetic  tapes.  Although 
the  BAMKI  format  was  able  to  record  all  the  sensor  data  in 
real  time,  it  caused  many  tape  read  problems  which  delayed 
processing  the  tape  at  the  Honeywell  635. 

While  the  BAMKI  format  offered  some  advantages,  it  also 
has  many  deficiencies.  The  time  required  to  unpack  the  45 
simultaneous  data  channels  made  the  BAMKI  format  unwieldy 
for  quick  access  and  analysis  of  the  data.  These  tapes 
contained  aborted  runs  which  should  have  been  deleted  but 
were  mixed  in  with  the  valid  runs.  Also  a  software  bug  in 
the  PDP-9  magnetic  tape  driver  resulted  in  a  high  rate  of 
parity  errors  when  we  tried  to  read  these  tapes  at  the 
Honeywell  635.  For  these  reasons,  the  seismic  data  was 
stripped  from  the  BAMKI  tapes,  edited,  and  formatted  more 
simply  on  other  tapes.  These  new  tapes  presented  the 
advantages  of  clean,  parity  error  free  data  and  a  simple 
format  which  made  the  data  easily  accessible. 
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SECTION  III 


SEGMENTATION 


In  order  to  analyze  "clean"  data,  i.e.,  data  which  is 
truly  characteristic  of  each  target  class,  a  criterion  was 
developed  to  cut  data  from  each  run  and  save  only  a  meaning¬ 
ful  portion  of  the  run.  This  segmentation  operation  is 
useful  because  (1)  it  presents  only  the  statistically  signi¬ 
ficant  data  to  the  decision  making  stage  of  logic,  thus 
promising  higher  recognition  rates,  ana  (2)  it  reduces  the 
amount  of  time  the  sensor  must  be  processing  data  for  deci¬ 
sions,  thus  reducing  power  requirements  and  extending  sensor 
life. 

Development  of  a  segmentation  procedure  requires  one 
major  step  in  common  with  feature  design:  extensive  visual 
study  of  the  vaveforms  on  hardcopy  or  graphic  displays.  In 
a  non-trlvlal  problem,  valid  features  can't  be  selected  and 
designed  until  the  engineers  have  a  very  thorough  knowledge 
of  the  signal  characteristics  of  each  class,  and  optimally 
u  thorough  understanding  of  the  physics  behind  these  charac¬ 
teristics.  This  in-depth  knowledge  of  the  data  should  allow 
the  design  of  a  reasonable  segmentation  algorithm  and 
criterions.  To  gain  this  required  Information,  signal 
waveforms  were  recreated  and  displayed  on  interactive 
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graphics  devices,  such  as  the  CDC  1704  Digigraphics  Display 
and  the  Tektronics  4002A  Graphics  Terminal.  Estimates  of 
the  energy  spectrum,  using  the  Fast  Fourier  Transform,  were 
computed  and  displayed.  Displays  and  hardcopies  of  these  wave¬ 
forms  and  their  power  spectrums  were  the  tools  which  enabled 
the  designers  to  view  and  analyze  the  behavior  of  each  class 
of  seismic  waveforms. 

Since  data  sementatlon  in  the  real  intrusion  detection 
system  will  probably  be  done  at  the  sensor,  simplicity  and 
efficiency  are  of  utmost  importance.  The  procedure  decided 
upon  begins  with  calculating  the  mean  value  of  the  entire  run, 
then  subtracting  that  mean  from  the  run  (realizable  in  the 
field  by  appropriate  capacitive  coupling  in  the  sensor's 
analog  output)  to  eliminate  any  DC  bias.  The  signal  is  then 
full-wave  rectified.  The  average  value  of  each  second  of  the 
rectified  signal  is  then  computed,  and  a  segment  of  valid 
data  is  defined  as  one  for  which  this  one-second  average  ex¬ 
ceeds  some  threshold  0  for  five  consecutive  seconds. 

Symbolically,  given  the  samples  f^  of  a  complete  in¬ 
trusion  run,  the  average  absolute  value,  S^,  will  be  calculated 
for  each  consecutive  one-second  window. 
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n  -  number  of  samples  per  one  second  window  and 
N  -  number  of  samples  in  complete  run 
Segment  and  save  the  data  in  the  five  second  interval  if  and 
only  if  Sk  is  greater  than  the  threshold  0  for  five  consecu¬ 
tive  one  second  windows. 


Obviously,  strong  signals,  from  either  large  sources  or 
intrusions  near  to  sensors,  will  result  in  a  greater  number 
of  five-second  segments.  This  is  desirable,  since  these 
stronger  signals  represent  a  better  s ignal-to-noise  ratio. 

The  specific  segment  lengths  and  thresholds  were  based 
on  observation  and  experimentation.  The  five-second  length 
precluded  the  acceptance  of  spurious  bursts  of  noise  or 
brief  signal  transients  as  good  data.  Also,  the  second-by- 
second  threshold  requirement  during  the  five  seconds  assured 
that  the  entire  segment  was  sufficiently  strong,  instead  of 
having  brief  but  significant  lapses  into  noise.  The  possi¬ 
bility  of  triggering  this  classifier  with  impulsive  noise, 
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such  as  explosions,  gun-fire,  etc.  is  not  likely  unless  the 
noise  were  highly  repetitive  and  sustained  over  a  five 
second  interval. 

The  threshold,  however,  luired  the  collection  of 
some  statistics.  The  objective  was  a  threshold  which  would 
overlook  as  much  noise  as  possible,  yet  which  would  locate 
as  much  valid  intrusion  data  as  possible.  We  selected  three 
representative  runs  from  each  data  class  (including  strictly 
noise  runs)  and  compiled  tables  of  the  total  time  segmented 
from  each  run  by  a  variety  of  thresholds.  We  then  observed, 
via  graphics,  the  five-second  segments  selected  by  thresholds 
of  20,  30,  and  40,  and  decided  that  6  *  20  afforded  the  best 
balance  between  noise  rejection  and  significant  data  segmen¬ 
tation. 
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SECTION  IV 


DATA  ANALYSIS 


Classifier  design  requires  Che  analysis  of  graphic  repre¬ 
sentations  of  digitized  waveforms  and  their  transforms  for  the 
purpose  of  hypothesizing  measurements  or  features  which  may 
aid  in  the  discrimination  of  target  classes.  This  is  one  of 
the  most  important  steps  in  the  waveform  classification  prob¬ 
lem  because  the  quality  of  the  selected  features  directly 
Influences  the  classifier's  performance. 

Before  data  analysis  can  begin,  the  researcher  must 
develop  or  have  access  to  a  system  which  will  display  his 
data  in  some  graphic  form.  At  the  start  of  ISCP's  Involvement 
in  the  sensor  program,  in-house  personnel  developed  a  waveform 
analysis  software  system  on  the  CDC-1700/Digigraphics  System. 
This  waveform  analysis  tool  enabled  a  user  to  randomly  access 
and  display  waveforms  or  waveform  segments  and  perform 
operations  on  the  data,  such  as  rectifying,  integrating, 
measuring  zero  crossings,  calculating  power  spectrums,  etc. 

The  major  deficiency  in  this  system  was  its  lack  of  a  hard 
copy  capability.  Unf or r unate ly ,  shortly  after  this  inter¬ 
active  graphics  software  system  was  operational,  the  CDC-1700/ 
Digigraphlcs  System  was  phased  out. 
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Effort  was  redirected  to  develop  a  similar  interactive 
capability  on  the  Honeywell  635,  using  a  remote  storage  tube 
terminal  as  the  graphics  console.  Soon  after  this  development 
started,  it  was  evident  that  the  development  time  required 
for  an  interactive  system  in  the  GECOS  III  multiprogramming 
environment  was  not  compatible  with  the  schedule  for  develop¬ 
ing  this  classifier. 

Pattern  Analysis  and  Recognition  Corporation  (PAR)  then 
made  their  sensor  analysis  factility  [3]  available  for  this 
effort . 

PAR's  sensor  analysis  facility  is  designed  to  analyze 
acoustic  data  and  is  built  around  a  NOVA  800  computer  with 
a  9-track  tape  unit,  card  reader,  128K  word  disk,  and  a 
Tektronix  4002A  display  with  hardcopy. 

Two  software  changes  had  to  be  made  before  we  could  use 
PAR's  facility.  First,  the  BAMKI  data  stripping  program 
was  modified  to  generate  its  output  on  the  Honeywell  635's 
only  available  nine-track  tape  drive  because  the  PAR  facili¬ 
ty's  only  tape  drive  was  a  nine-track  unit.  Second,  PAR 
modified  their  system  by  adding  a  new  input  routine  to  it, 
which  could  read  the  nine-track  version  of  our  simplified 
data  format. 

Once  we  were  able  to  process  seismic  data  from  a  BAMKI 
tape,  using  PAR's  facility,  a  production  procedure  was  set 
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u  .  Three  copies  of  Che  BAMKI  Capes  could  be  lefC  Co  be 
processed  aC  nlghC  by  Che  635.  Since  Che  Honeywell  635  has 
only  one  available  nine-Crack  Cape  drive,  only  one  copy  of 
Che  BAMKI  scripping  program  was  able  Co  execuCe  aC  a  Cime. 
Ofcen  Che  BAMKI  programs  were  delayed  from  execudng  for 
long  periods  of  Cime  because  Che  635's  only  nine-Crack 
Cape  drive  was  previously  assigned  Co  long  baCch  jobs. 

Time  domain  waveforms  and  power  specCrums  for  at 
Least  one  complete  run  of  each  intrusion  class 
variation  were  displayed  and  hardcopied  aC  Che  sensor 
analysis  facility.  For  ground  vehicles,  these  class  varia¬ 
tions  were  Che  different  vehicle  velocities  recorded  along 
each  of  these  paths.  For  aircraft,  these  class  variations 
were  the  different  altitudes  of  the  flyovers  and  the  differ¬ 
ent  velocities  recorded  for  each  altitude.  The  class  varia¬ 
tions  for  humans  were  in  the  number  of  intruders,  path  of 
intrusion,  and  velocity  (feet  per  second).  Examples  of 
time  waveform  and  power  spectrum  hardcopies  for  each  class 
are  shown  in  Figures  1  thru  10. 

After  the  hardcopies  of  the  selected  intrusion  runs 
were  generated,  they  were  added  to  the  hardcopy  library. 

This  library  consisted  of  two  note  books,  one  containing 
the  time  domain  waveforms  and  the  other  containing  power 
spectrum  plots. 
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A  team  of  ISCP  engineers  analyzed  the  data  plots  con¬ 
tained  in  the  hardcopy  library.  By  comparing  the  various 
plots  of  the  intrusion  classes  looking  for  within -class 
similarities  and  between-class  differences,  this  team 
compiled  a  list  of  48  potential  features.  One  very  useful 
data  plot  used  in  our  analysis  was  the  40-line  power 
spectrum  in  which  each  line  represents  the  power  spectrums 
of  consecutive  one-second  windows.  The  representation  of 
the  data  in  this  format  enables  the  analyst  to  view  and 
compare  changes  in  the  power  spectrum  throughout  the 
duration  of  the  intrusion.  Examples  of  the  40-line  power 
spectrum  are  shown  in  Figures  11  and  12.  Figure  11  shows 
the  40-line  power  spectrum  plot  of  a  C-131  aircraft  on  a 
radial  path  over  the  test  site.  The  frequency  shift  in 
the  main  peaks  of  the  power  spectrum  indicates  when  the 
C-131  responses  went  through  a  doppler  shift.  A  similar 
data  display  for  the  UU-1F  helicopter,  shown  in  Figure  12, 
indicates  that  the  doppler  shift  is  not  as  pronounced  for 
the  helicopter  as  it  was  for  the  C-131. 


CHA 


FIGURE  1 
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SAMPLE  WAVEFORM  OF  ONE  MAN  RUNNING 


SAMPLE  POWER  SPECTRUM  OF  ONE  MAN  RUNNING 


ale:cm  input:  :sex9.  u  lee  arc  irhphcpa  2«i  vert 


FIGURE  3 
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SAMPLE  WAVEFORM  OF  AN  ARMOURED  PERSONNEL  CARRIER 


FIGURE  4 
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SAMPLE  POWER  SPECTRUM  OF  AN  ARMOURED  PERSONNEL  CARRIER 


SAMPLE  WAVEFORM  OF  AH  M-109  TRUCK 
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SAMPLE  POWER  SPECTRUM  OF  AN  M-109  TRUCK 


*lc:lu  input ! H  jseis.  uh-if  SERP  75MPH  5®FT 


HGURE  7 
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SAMPLE  WAVEFORM  OF  A  UH-IF  HELICOPTER 
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SAMPLE  POWER  SPECTRUM  OF  A  UH-1F  HELICOPTER 


FIGURE  9 


20 


- 

- 

- 

- 

m 

- 

— 

- 

o 

■t 

wmm 

mm 

w 

— 

— 

— 

— 

MB 

►» 

• 

a. 

►- 

X 

__ 

m 

M 

• 

♦ 

1  ♦ 

r» 

mm 

♦ 

-u 

“  u 

Ul 

“'•B 

UJ 

on 

♦ 

Ul 

w. 

rv 

CO 

• 

Oi 

*4 

• 

n 

>.  <x 

—  n 

— •* 

m 

X 

X 

X 

X 

< 

« 

« 

c 

X 

X 

"  3C 

X 

(ft 

♦ 

~m» 

—  M> 

—  ♦ 

HO 

♦* 

v« 

ui 

Ul 

Ul 

Ul 

W 

_  PI 

_  a> 

r» 

CD 

UD 

rs 

CD 

m 

«4 

fN 

u 

CO 

♦ 

« 

O 

-« 

— 

— <* 

u 

■  1 

I- 

*- 

h- 

B- 

_o 

—  o 

_  O 

— O 

►- 

►- 

I- 

B- 

3 

• 

-v 

-  CD 

~  « 

— XX 

m 

«• 

CD 

o> 

•* 

m* 

«■ 

M  «  • 

—  •• 

B^»* 

k. 

'■*  ‘  ifl 

«n 

«* 

•• 

E 

u 

pi 

♦ 

“  Ul 

“Ul 

~u 

C 

K 

i  »- 

M  < 

-  h 

1  • 

i 

_♦ 

_  ■* 

__  ♦ 

X 

X 

X 

X 

z 

« 

« 

« 

<r 

« 

X  « 

-x  — ' 

—  X  « 

—  X 

“X 

•• 

»•  «_> 

o 

u» 

Ul 

Cl 

K 

3 

_ 

mmm 

MB 

wm 

a. 

a.  in 

r\ 

rv  4 

IN  < 

IN 

X 

jx 

** 

Hi  Z 

z 

X 

y 

X 

3 

3 

“  3 

■“O 

a  - 

SLJpH 

“ _ 

X 

W-  - 

>M 

♦ 

liT 

♦ 

♦ 

mt 

w4 

n 

M  -i 

“  <M 

•• 

M» 

K» 

ID 

3 

v 

—  x 

— *x 

U 

i  u 

i  <-> 

c_» 

Ul 

* 

►- 

ir _ 

Lr_ 

B- 

1- 

FIGURE  ] 


21 


SAMPLE  POWER  SPECTRUM  OF  A  C-131  AIRCRAFT 


*lc:lp  input:  :se!9.  u  lee  uh-if  r/o  ii»hph  >tt  veit 


SAMPLE  40  LINE  POWER  SPECTRUM  OF  A  UH-1F  HELICOPTER 


SECTION  V 


FEATURE  DEFINITION 


Several  man-months  of  studying  seismic  waveforms  on 
graphics  and  hardcopy  displays  produced  a  set  of  24  features 
for  consideration.  These  features  were  extracted  from  one- 
second  (actually  1024  samples,  or  1.024  seconds)  segments  of 
data  and  from  five-second  segments  (5120  samples,  5.120 
seconds)  of  data.  The  two  different  lengths  were  selected 
to  evaluate  the  effect  of  segment  length  on  classification 
success.  The  FFT  routine  required  the  number  of  samples  used 
to  ke  a  power  of  two.  Therefore,  the  one-second  window  size 
waachosen  to  contain  1024  samples  instead  of  1000.  These  24 
features  taken  over  the  two  segment  lengths  produced  a  total 
of  48  features. 

The  rationale  and  definitions  of  the  features  are  given 
below,  in  the  order  in  which  they  appear  in  the  vector  data. 
That  is,  component  1  of  the  vector  is  the  average  R  for  one 
second.  Any  DDC  offset  present  was  subtracted  before  all 
processing. 

NO.  DEFINITION  AND  RATIONALE 

1.  R  for  one  second  (Average  R  for  each  half- 

second  of  a  contiguous  one-second  segment): 

R  is  defined  as  the  ratio  of  the  maximum 
absolute  signal  amplitude  during  a  half- 
second  interval  to  the  average  absolute 


NO.  DEFINITION  AND  RATIONALE  (continued) 

amplitude  for  that  half  second.  This 
feature  appeared  a  likely  candidate  for 
distinguishing  the  class  of  humans,  whether 
one  or  several,  walking  or  running.  The 
Impact  of  a  heel  gives  a  sharp,  strong 
spike  in  the  signal,  which  decays  to  noise 
level  considerably  before  the  next  heel 
Impact.  This  effect  produces  a  significantly 
higher  value  of  R  than  does  a  vehicle  or 
aircraft,  since  the  latter  usually  produces 
high  spikes  only  when  the  signal  is  strong 
enough  to  produce  a  high  RMS  value. 

2.  R*  for  one  second  (average  Ra  for  the  two 
half  seconds):  Although  any  information 
contained  in  R  (e.g.,  a  threshold)  will  map 
uniquely  into  R2,  squaring  R  made  the 
distinction  between  humans  and  all  other 
classes  more  obvious  to  the  operator,  and 
eventually  to  OLPARS. 

3.  R  for  five  seconds  (average  R  for  10  contigu¬ 
ous  half-second  segments). 

4.  R2  for  5  seconds. 
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DEFINITION  AND  RATIONALE  (continued) 


NO^ 

5.  Harmonic  spacing  for  one  second  (the  most 
frequently  occurring  pairwise  spacing  between 
the  six  largest  peaks  in  the  power  spectrum 
above  40  Hz):  This  feature  was  suggested  by 
the  evenly  spaced  harmonics  that  were  evident 
in  aircraft  waveforms. 

6.  Harmonic  occurrences  in  one  second  (the  number 
of  times  that  the  spacing  of  component  5  above 
occurred  in  the  one  second  spectrum):  The  C-131 
usually  evidenced  harmonics  at  20-Hz  inter¬ 
vals,  the  UH-1F  helicopter  at  12.5  Hz  inter¬ 
vals. 

7.  The  ratio  of  the  energy  between  1  and  20  Hz 
to  that  between  21  and  40  Hz  for  one  second 
(hereafter  symbolized  by:  Ei-20^E2l-40^  5 
Although  the  jeep  and  APC  spectra  overlapped 
considerably,  the  jeep  spectrum  did  extend 
somewhat  lower  than  the  APC  spectrum,  which 
usually  dropped  off  below  20  Hz. 

8.  E4 1-60 ^E2 1-40  ^or  one  second!  An  attempt  to 
separate  trucks  from  APCs. 
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NO.  DEFINITION  AND  RATIONALE  (continued) 

9.  The  number  of  points  in  the  one  second  power 

spectrum  which  are  below  25 %  of  the  maximum: 

A  coarse  amplitude  histogram  was  taken  to 
estimate  the  value  of  amplitude  information, 

10.  The  number  of  points  in  the  one  second  power 
spectrum  which  are  2125%  and  <502  of  the 
maximum . 

11.  Number  as  in  9  and  10  of  points  >^50%  and  <75%. 

12.  Number  of  spectral  points  _> 75%  of  maximum. 

13.  The  ratio  of  the  energy  above  100  Hz  to  that 
below  in  the  one-sec.  power  spectrum ( i . e  . 

E101-  'ill/El-100):  Aircraft  tend  to  produce 
more  energy  above  100  Hz  than  do  any  other 

classes . 

14.  El-5^El-60  for  the  first  second:  To  estimate 
the  value  of  spectral  distribution  information. 


15. 

E6-10^El-60 »  first 

second . 

16. 

Ell-15^El-60 ’  first 

second . 

17. 

E16-20/,E1”60 »  first 

second. 

18. 

E21-25^El-60 »  first 

second . 

19. 

e26-30^e1-60 »  first 

second . 

20. 

• 
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No 


DEFINITION  AND  RATIONALE  (continued) 

25.  *“56-60^1-60  ’  ^*r8t  second. 

26.  The  number  of  peaks  in  the  first  one  second 
power  spectrum  which  exceed  10%  of  the  maximum 
peak. 

27.  Harmonic  spacing  for  the  5-second  ensemble 
average  of  5  successive  one-second  power 
spectra  (see  No.  5). 

28.  The  number  of  times  the  spacing  of  No.  27 
occurs  in  the  5-second  ensemble  average 
(see  no .  6">  . 

29-A8.  Components  /9-A8  are  extracted  just  as 

components  7-26,  in  that  sequence,  except 
that  the  spectra  observed  consist  of  an 
ensemble  average  of  5  consecutive  one- 
second  power  spectra  instead  of  a  single 
one-second  spectrum. 
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SECTION  VI 
FEATURE  EXTRACTION 


Once  Che  Use  of  48  features  had  been  compiled*  a  batch 
program  was  written  for  the  Honeywell  635  to  extract  these 
features  from  the  unpacked  data  tapes.  This  batch  program 
consists  of  the  main  "control"  program*  a  Fast  Fourier  Trans¬ 
form  (FFT)  subroutine*  a  double  up  algorithm,  and  a  number  of 
feature  extraction  subroutines. 

Data  cards  direct  the  control  program  as  to  which  and 
how  many  runs  and  channels  are  to  be  processed.  The  control 
program  then  monitors  the  specified  input  data  tape  channels 
and  calls  the  feature  extraction  subroutines  when  the  moni¬ 
tored  data  channel  satisfies  the  automatic  segmentation 
criteria.  The  first  subroutines  called  extract  features 
from  the  time  domain  waveform  segment  in  the  data  array  X. 
After  these  features  are  extracted,  the  double  up  algorithm 
is  called  three  times  to  help  calculate  five  consecutive 
one-second  power  spectra  of  the  data  in  X.  Since  the 
waveform  segments  do  not  have  an  imaginary  component,  the 
double  up  algorithm  enables  the  FFT  to  calculate  two  power 
spectra  with  one  call  to  subroutine  DOUBLE.  After  the  power 
spectra  are  returned,  their  ensemble  average  is  computed. 
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The  X  array  will  then  contain  two  power  spectra  of  interest, 
the  spectrum  for  the  1st  second  and  the  average  spectrum  for 
the  five  seconds  of  data.  The  subroutines  which  extract 
features  from  these  two  power  spectra  are  then  called 
After  all  the  feature  extraction  subroutines  are  called,  a 
labeled  feature  vector  is  punched  out  on  cards. 

After  completion  of  the  feature  extraction  program,  the 
output  data  deck  is  taken  to  the  CDC-1604B  computer  where  the 
feature  vectors  are  transferred  to  an  OLPARS  compatible  tape. 

Listings  of  the  feature  extraction  program  and  its 
subroutines  are  provided  in  Appendix  A. 
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SECTION  VII 


FEATURE  EVALUATION 

The  first  step  in  using  the  On-Line  Pattern  Analysis 
and  Recognition  System  (OLPARS)  is  to  evaluate  the  discrimi 
natory  quality  of  the  extracted  features.  This  enables  us 
to  use  fewer  measurements  to  achieve  a  satisfactory  classi¬ 
fier  design.  The  OLPARS  provides  two  suboptimal  methods 
for  ranking  the  discriminatory  power  of  the  extracted 
features.  Each  of  these  methods  provides  three  types  of 
rankings.  The  first  type  uses  a  s ignif i cance  measure  of  a 
particular  feature,  Xp ,  for  discriminating  class  i  from 
class  j  and  is  designated  by  (Xp) .  The  second  type  of 
ranking  uses  a  significap.ce  measure  of  Xp  for  discriminating 
class  i  from  all  other  classes  and  is  designated  Mi(Xp). 

The  last  type  uses  a  measure  of  the  overall  significance 
of  Xp  for  discriminating  all  classes  and  is  designated 
M(Xp)  . 

The  first  method  in  OLPARS  for  ranking  features  is 
the  discriminant  measure,  which  is  particularly  useful  when 
the  class  conditional  probability  distributions  are  unimoda 
These  discriminant  measures,  using  feature  Xp,  are  defined 
as  f  o 1 lows : 
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Mij<Xp) 


y(n  "y  (  j  )“I2 

*  p  *  p 

(..-dCs:1]*  + 

K 

Mi(Xp)  “  X  Mij<Xp) 

tC  K  K 

«  (X  )  ■  S  «1<X  )  -  X  Z  MU(Kp> 
i«l  i»l  j*i 

where  Xp^  *  the  estimated  mean  of  class  j  along 

measurement  Y  • 

P 

A  (  1  ) 

Op  ■  the  estimated  standard  deviation  of 
class  J  along  measurement  xp* 

N  ^  =  the  number  of  samples  from  class  j. 

The  other  OLPARS  feature  evaluation  method  is  the 
probability  of  confusion  measure.  It  is  valid  for  any 
probability  distribution  since  it  essentially  measures 
the  overlap  of  the  class  conditional  probabilities. 

Since  the  functional  forms  of  the  class  conditional 
probabilities  are  not  known,  OLPARS  estimates  the  marginal 
class  distributions  using  the  sample  data.  The  range  for 
feature  Xp  is  divided  into  cells  of  width  A.  The  probabili¬ 
ty  that  a  sample  from  class  1  will  occupy  the  cell  along 

the  range  of  feature  Xp  Is  f.  i  v  e  n  by: 
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where  N  -  tiie  number  of  cells  alone  measurement  Y„ 

p  P 

and  K  =  die  number  of  classes. 

Tlie  ranking  of  extracted  features  based  on  these  eval¬ 
uation  techniques  provides  the  information  required  to 
rationally  choose  initial  subsets  of  the  48  features  for 
logic  design.  Logic  design  is  an  iterative  process  in  which 
many  designs,  based  on  modified  versions  of  t he  initial 
feature  subsets,  are  generated  and  tested.  Features  which 
appear  to  discriminate  between  the  more  troublesome  classes 
are  added,  while  superfluous  features  which  rank  high  for 
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the  same  easily  discriminated  classes  are  eliminated. 

For  this  five-class  problem,  the  top  fifteen  features, 
rank-ordered  by  the  probability  of  confusion  measure,  ^ij(X  ), 
are  shown  in  Appendix  B. 
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SECTION  VIII 
CLASSIFICATION  LOGIC 

The  classifiers  designed  by  IS CP  for  this  pattern  recog¬ 
nition  problem  consist  solely  of  sets  of  linear  discriminants 
for  ease  of  hardware  implementation.  The  logic  for  these 
classifiers  is  based  on  the  pairwise  Fisher  Linear  Discrimi¬ 
nant  Technique.  For  each  pair  of  classes  i  and  j,  a  unit 
vector  dij  is  computed  such  that  projections  of  the  data 
onto  d^j  maximize  the  ratio  of  the  between-class  scatter  to 
the  within-class  scatter.  The  direction  d^j  which  maximizes 
this  ratio  is  given  by  Reference  [5]. 

its  m  °  wij 

where  WiJ  -  (Ni  -  1)  Cj  +  (Nj  -  1)  Cj 

C^j  -  Estimated  covariance  matrix  for  class  i 

Aij  “  Hi  -  Hj 

Ui  -  Estimated  mean  vector  of  class  i 

N^  ■  Number  of  vectors  in  class  i 
and  a  is  a  normalizing  constant  so  that  |d|  =  1 
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OLPARS  computes  d^j  and  an  Initial  threshold,  O^j ,  to 
distinguish  between  all  pairs  of  classes.  These  thresholds 
may  be  adjusted,  if  necessary,  to  obtain  optimal  discrimina¬ 
tion  along  each  d^j. 

For  example,  the  inner  product  of  an  unknown  input 
feature  vector,  3,  is  taken  with  the  discriminant  dAH  for  the 
pair  of  APCs  and  helicopters,  compared  with  the  threshold 
®AH  ^or  Pa^r  of  APCs  and  helicopters. 

increment  the  counter 
for  the  APC  class. 

increment  the  counter 
for  the  helicopter  class. 

increment  the  counter 

4 

for  the  class  with  the 
larger  number  of  samples 
in  the  design  set. 

where  K  =  Size  of  the  feature  space. 

After  all  the  pairwise  decisions  are  made,  a  binary  vote 
is  cast  by  each  comparator  and  the  final  decision  is  deter¬ 
mined  by  the  class  counter  that  received  the  most  votes.  In 
case  of  ties,  the  decision  is  given  to  the  class  involved  in 
the  tie  which  has  the  highest  a  priori  probability.  The  re¬ 
sultant  classification  scheme  is  diagrammed  in  Figure  13. 


If  <^AH*X>  "X  Xi  *AII  >0Ail 
i-1 


If  <  -All  *  “  X  Xi  -All  <0AH 

i-1 


If  <iiAlI»X>  “  Z  Xt  *2 AH  -  0AU 
i-1 
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PAIRWISE  FISHER  LOGIC 
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Figure  13, 
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h  ■  UNKNOWN 
INPUT 

INVERTER 

d  >  LINEAR 
DISCRIMINANT 

e>  THRESHOLD 


H'  HELICOPTER 
A  :  APC 

T  .  HEAVY  TRUCK 
M<  HUMAN 
C  C- 131 


SECTION  IX 


RECOGNITION  RATES 


Four  classifiers  were  designed  using  16,  2  2,  33,  and 
44  features.  With  the  exception  that  feature  46  is  missing 
from  the  44  feature  design,  each  design  is  based  on  a  subset 
of  the  features  used  in  the  higher  order  designs.  Figure  14 
lists  the  features  used  for  each  design. 

These  classifiers  were  evaluated  with  the  design  data 
set  and  an  independent  test  data  set.  The  design  set  con¬ 
sisted  of  715  vectors:  133  vectors  from  the  class  of 
helicopters;  176  for  APCs;  135  for  heavy  trucks;  136  for 
humans;  and  135  represented  C-131s.  The  independent  test 
set  was  comprised  of  607  vectors:  159  for  lie  1  icopters  ; 

171  for  APCs;  16  for  heavy  trucks;  129  for  humans;  and  132 
for  C-131s.  The  design  and  test  confusion  matrices  from 
the  resulting  evaluation  of  each  classifier  are  shown  in 


Figures  15  thru  18. 
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CONFUSION  »IATRIX  FOR  THE  DESIGN  SET 
USING  22  FEATURES 

ASSIGNED  CLASS 

HAT  M  C 

II  128  4  5  4  18 

A  3  144  20  4  0 

TRUE 

CLASS  T  1  3  12  0  0 

MO  5  0  123  1 

C  11  15  0  115 

Probability  of  Correct  Classification  *  .860 
CONFUSION  MATRIX  FOR  THE  TEST  SET 
USING  22  FEATURES 

LEGEND:  II  -  HELICOPTER  M  -  HUMAN 

A  -  APC  C  -  C-131 

T  -  HEAVY  TRUCK 

FIGURE  17 
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ASSIGNED  CLASS 


TRUE 

CLASS 


H 

A 

T 

M 

C 


II 

120 

3 

1 

2 

1 


A 

2 

170 

4 

1 

n 


T 

4 

2 

130 

0 

2 


M 

4 

1 

0 

120 

1 


Probability  of  Correct  Classification  *  .OS] 
CONFUSION  MATRIX  FOR  THE  DESIGN  SET 


C 

3 
0 
0 

4 

131 


USING  16  FEATURES 


ASSIGNED  CLASS 


II 

A 

T 

M 

c 

II 

120 

2 

5 

4 

19 

A 

6 

142 

19 

4 

0 

TRUE 

CLASS 

T 

1 

3 

12 

0 

0 

M 

1 

6 

n 

121 

1 

C 

11 

1 

6 

n 

114 

Probability  of  Correct  Classification  *  .853 
CONFUSION  MATRIX  FOR  THE  TEST  SET 


USING  16  FEATURES 


LEGEND:  II  -  IIEIICOPTER  M  -  HUMAN 

A  -  ARC  C  -  C-131 

T  -  HEAVY  TRUCK 


FIGURE  18 
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Note  that  the  probability  of  correct  classification 
using  the  design  set  increases  monotonically  as  the  number 
of  features  is  increased;  while  using  the  test  set,  it 
reaches  its  maximum  value  for  22  features  and  drops  off. 
Foley  [2]  points  out  a  statistical  trap  involved  with  using 
the  probability  of  error  on  the  design  as  a  measure  of  the 
true  performance  of  the  system  when  the  ratio  of  the  sample 
size  to  feature  size  is  small.  The  probability  of  error 
or  correct  classification  on  the  test  set  is  the  better 
measure  and  indicates  that  the  probability  of  correct 
classification  of  0.86  for  the  22-feature  design  is  the 
best  performance. 

The  essential  difference  between  the  22-feature  design 
and  the  16-feature  design  is  that  the  former  is  based  on 
measurements  of  both  one  and  five  second  data  segments, 
whereas,  the  latter  is  based  solely  on  measurements  of 
five  second  segments.  Therefore,  the  exclusion  of  one 
second  measurements  not  only  reduces  the  dimensionality 
of  the  decision  logic  from  22  to  16,  but  it  also  simplifies 
feature  extract'. ion  by  elimination  of  the  one  second  phase 
of  the  extraction  process.  The  trade-off  between  the 
slightly  superior  performance  of  the  22-feature  design 
(Pcc  on  the  test  set  of  0.86)  versus  i'ne  simplicity  and 
ease  of  implementation  of  the  16-feature  design  clearly 


ranks  the  16  feature  design  ahead  of  the  22-feature  design. 

The  linear  discriminants  for  the  16- feature  logic 
design  are  given  in  Appendix  C.  These  linaar  discriminants 
contain  the  ten  weight  vectors  with  their  respective  thresh 
olds  and,  together,  with  the  block  diagram  in  Figure  13, 
completely  define  the  pairwise  Fisher  logic  for  the  16- 
feature  design. 
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SECTION  X 


DISCUSSION  AND  RECOMMENDATIONS 


In  many  waveform  classification  problems,  the  classifier 
is  designed  entirely  from  features  extracted  from  a  waveform 
data  base  which  is  representative  of  each  class.  To  insure 
that  the  data  base  is  representative,  data  must  be  collected 
for  a  sufficient  number  of  runs  of  each  intrusion  variation 
(speed,  path,  direction,  etc.),  so  that  the  data  for  each 
variation  is  truly  representative.  Emphasis  must  be  placed 
on  the  generality  of  a  data  base  and  not  on  its  size  alone. 
For  example,  to  design  a  classifier  which  will  detect  most 
trucks,  the  data  base  must  contain  data  from  a  wide  variety 
of  different  type  trucks,  all  operating  at  a  number  of 
sampled  speeds  and  carrying  loads  which  vary  from  maximum 
capacLty  to  empty.  The  wider  the  variety  of  vehicles  within 
a  class,  the  more  complex  the  data  analysis  becomes  because 
all  the  variations  and  the  possible  combinations  of  the 
variations  have  to  be  sampled  and  analyzed  for  each  variety 
of  vehicles. 

The  mair  goal  of  this  effort  is  the  development  of  a 
seismic  classifier  which  w.r».l  satisfy  an  acceptable  error 
criteria  regardless  of  where  the  classifier  is  located. 

The  seismometer  response  to  an  intrusion  target  is  the  signal 
which  results  from  the  convolution  of  the  target  signature 
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with  the  impulse  response  of  the  earth.  The  earth's  impulse 
responses  differ  greatly  throughout  the  world  because  they  :.re 
determined  by  the  geology  and  terrain  of  the  locale.  These 
impulse  responses  or  transfer  characteristics  are  not  neces¬ 
sarily  constant  between  two  points  in  a  locality  but  may  change 
periodically  as  a  function  of  the  seasonal  variations.  Sea¬ 
sonal  variations  in  the  transfer  characteristics  are  caused 
by  changes  in  either  the  height  of  the  water  table,  the  depth 
of  the  frost  line,  or  amount  of  snow  cover.  Therefore,  it 
may  not  be  realistic  to  believe  that  one  seismic  logic  design 
might  perform  satisfactorily  in  any  location  or  for  an  extended 
period  of  time  in  a  locality  which  has  extreme  seasonal 
variations . 

The  key  to  the  general  seismic  classifier  design  problem 
is  the  discovery  of  waveform  features  which  contain  adequate 
discriminatory  information  and  are  invariant  to  geographical 
location  and  seasonal  variation.  This  report  does  not  shed 
any  light  on  the  theory  that  such  features  exist,  but  it 
does  exercise  the  procedures  required  to  determine  their 
existence.  The  importance  of  the  data  collection  and  data 
analysis  procedures  cannot  be  over  emphasized.  Before  data 
analysis  can  yield  meaningful  results,  the  data  base  must 
contain  data  collected  from  selected  test  sites  which  are 
representative  of  the  various  types  of  terrain,  geology,  and 
seasonal  variations  found  throughout  the  world. 
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appendix  a 


BATCH  PROGRAM  LISTINGS 
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0 4 -0 4 " 7  3  17,594  PROGRAM  TO  STRIP  SPECIFIED  DATA  CHANNELS  FROM  A  BAMK1  TAPE 

_  AND  WRITF  A  Ngh  TAPE  WITH  A_  SJ  MPL  I  P  I.ED  FORMAT, 

C  PROGRAM  TQ  STRIP  SPECIFIED  DATA  CHANNELS  FROM  A  RAMK!  TAPE 

0  AND  WRITE  A  NEW  TAPE  WITH  A  SIMPLIFIED  FORMAT, 

CQKMCN  I  OAT  <  3i>U »  V )  #  JCOUNT  <  45  ) 

COMMON  /JAZZ/lAIjbSB)  , I  EOF ,  I  BORT , I  PART Y 

_ UMAQN/JIM/LAL326) _ _ _  _ _ 

C0MMCN/AVbBGP/IAVG(45) ,lC0UNT(45)#  jjk.NTIMES, jPRiNT 
DIMENSION  KCOyNT (45) 

IPARTYsG 
IBCRT  *C 
I  JK*o 

R£Ar)  UOl.NJM,  (KcOUNTtjJjjsl.NUM) 

1001  f  Of-  MAT  <  40  I  2  ) 

READ  1003.NREC.NTIMES 

1003  F  OF  M  A  T  (  2  l  5  ) 

READ  100i.IPRlNT 
LU  4 L  0  I  cl i NUM 

_ imuKCflUNim _ _ _ 

wR1TF(6,1006  )  I  *  I ND 

1006  fORMaT(10X.7hRfCORd  , I  2 • 1 7m  IS  fROM  cHaNN^L  .1?) 

4  0  G  L  Cg’Jn  T  C  ImD  )  *  T 

DO  3  J«l,45 
3  JCCJNT(J>sn 

_ - _ LO  _lilD5  M*liNRK _ _  _ _ __  _ 

CALL  SPAWNIN.NUM.KCOUNT ) 

I  F (  It  OF  . EQ . 1 >  GO  TO  1004 
1F<  IFoPT.NE.O)  GO  To  1005 
I F  c  IjK.NF.O)  GO  TO  1005 
LO  1C00  Jsl.MUM 

_ _ 

LO  510  1  *  7  »  3 ?6 
500  LA(l)»lD*T(I-o.J* 
c ah.  ppoc 

1030  CONTINUE 
1005  CONTINUE 

1004  irtUK.EQ.P)  GO  TO  1007 _ 

I F (  IFOPT.NF.O)  GO  TO  1007 

LO  65  1*1,  N'JM 
INC*KC0UNT( I ) 

KEI'  D*  JOOUNT  (  I  NO ) 

LO  6 t  jxKgND • 320 

66  _J  D  A  T  (  J ,  L)  c  0 _  ________ 

L  A  (  4  )  *  I  A  VG  (  1  ) 

LO  67  J=7 , 326 

67  L  A  <  J ) *  I D*  r  f  J-6 »  I  ) 

65  CALL  PROC 

1 0 0 7  CALL  wFF 

_ CALu  WFF  _  _  _  _  _  _ 

call'wef 

STOP 

END 
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01  04-0-1-73 


1 1. .  ^  1 


I 


SUb^CJTjf  6  5>P  A  (  jP  A?m  ,  NUM  »  ft  CO'*NT  ) 

C 

C(JM>lC.\/Jlrt/t.A(  32*0 

COMf-N  /JA7^/IA(ld‘>t>)  .  IE0F.  Ifci^KT,  tPAATV 
COMMON  I  P  A  T  (  3  2  y  ,  7),  JCPU\T<45> 

C0MM0\/AVfc'5(,F/I  A  VS (45  ' » LCO'jNT  ( 46  ) »  1  JK  ,  *'T  !  Mt  S»  /  p‘  R  I  <4  T 
blMEf-SIOf  I  On  (  4  ) ,  JOL'  <  J?  >  •  JP  <  14  > 

L  If;Ef>  S  IQI"  xCnUMT  (  43)  ,  m  AS1'  ( 3 )  .  KPtvST  (  3  ’ 
i  1  KENS  10*  FAVoM*) 

LATA  I«ALF /2621 44/ 

L  AI  A  (  JP(  l ) ,  Isl,1  4>/l  ,?1 ,1  ,l?,?2»3n.i  .  1  ?.1.12,'<1 ,39,40.45/ 

LATA  K*’ /P/77777  7  70900/ 

[ATA(^ASK(l),I  =  l,3)/O7  7770nC)0C000,n0''0r77770OCl,n0n0O0r00  /777/ 
LATA  ( Ky A bT <1), I ri, 3)71677/216. 40 9ft, 1/ 

LATA  MASKA,HASKb/07/77770CC'ODC,GOOPln'l'777777/ 

40  F  OriMAT  (6(2/,  1 1?)  ) 

LALL  tf^ll KP 
IP(  n-pr.f’0.1)  P£TU»v 
I r ( JF  ASM  .NF  .1 )  6P  TO  2 
LA(3)sANP(  I  A  (  2  >  ,N‘ASKA)/JHALF 
L  A  0  4  C  sAM!(  I  A  (  2)  ,  VASftP  ) 

I S  A  ^  F i * TlrtFS 
NT IMFSslOOO/ TbAMF 
Ff  COF  DrO^_ 

?  10LU  )sAf  J(  14(1),  MASK  A)/ I  HALF 

I0i)<2  )  a  ANL)  (  I  A  (  J_  )  ,  MaKfc  ) 

10j(3>  =  AMU(  IA(?)  ,MASKA)/!H/LF 
l oy(  4 )  =  a;>'u<  i  a ( ? > , k a sk b ) 

IF(LA04C.NP .100(1))  00  TO  7 
I  F  (  L  A  (  3  J  ,_b  0  ,  I  0 OJ  3 )  )  00  TO  5 
i  F  (  HOST)  6, 19, 10 
5  IF  (  JFGr<T)6«b«5PCr 
ft  CALL  hSF 
CALL  FPp 
IfeORTsl 

GO  TC  5Q00  _  _  _ 

7  'iF  (  LA  (  7  >  TfcT,  I  on  (  3  >  )  00  TP  11 
I  F  (  IF  0»T  Ifc.jQ.TO 
11  I F ( I) OR T) 6,1?, 1C 
19  I  F  <  Ijk.Ec,  n)  on  Ty  9 
10  6b  1  =  1  .N'UM 
_|\LsK(;PU^'T  (  I  ) 

KEND«JC0(JNT(IND) 
yO  6ft  J  =  K£M[),3?0 
ft6  I  0 AT  <  J  ,  I  )  =  0 

L  A  (  4  )  s  I  A  V  0  (  I  ) 

LO  67  J * 7 , 3 2 ft 
67  L*< J)=IPAT( j-6, I) 

65,  CALL'  PPdfi  ' 

00  TC  9 
1.2  CALL  BSF 
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11  04-03-73  U.S15 


_? _ CALL  M£E _ 

10  1 0ORT  *  0 

RECORD *0. 

1JKIL 

9  LAI  1 ) * IOC( 1 ) 

LA  ( 2 ) * !  S  AMP 

_ LA ( 3  ?  » I  On  <  3  ) _ 

L  AO4C  * ! on ( 4 ) 

t.A(5)sANn<  !A(226)  »MASKP) 
LA(6)*ANn( 1  A (229) ,  MSK  A  )  / 1 H  ALF -4  0 
l F <  L A  <  6 ) , GP 1 0 )  GO  TO  3 
L  A(5)«LA(*>)-1 

_  _  LA(6)«1000»LA(») _ _ 

3  I F <  I  JK.GT.O)  GO  TO  444 
LO  50  J«1 » NUM 
InD*K  CCUmT ( j ) 

5Q  JCOJKT ( I  NO ) >1 
444  I C*  2 


LO  8 C 0  MMsi.io 
00  400  K*1 » 9 »  2 
«lBE*  JP  (  K ) 
jPt* JP ( K*1 ) 

KtO 

_ L£«_L£*_1 _ 

DO  400  J*J0E.JPE 

I F ( M , L T . 4 )  GO  TO  401 
1C* I C ♦ 1 
Ml 

401  IF ( JCOUNT ( J) , EO  «  0  )  GO  TO  4 0 0 
1 1* AND ( I  A ( ! C ) 1 MASK I  M )) 
ll«!l/KONST(M) 
lFdl.GE.2n48)  1 1  sOH  (  1 1  #  KM  ) 
Jl« JCOUNT ( J ) 

H«LC01JNT<  J) 

1 DaT ( J1 >  LI ) ■ 1 1 _ 

JCOUKT  <  J) * JCOUnT  <  J) M 
400  CONTINUE 

IF(MK-2)600, 700.800 
600  ^Bfc* JP ( 11 ) 

JPE«JP(12) 

JjO  _TO  4fi| 

700  JBE*JP(1J) 
jPfc«jP(14) 

402  *«0 
IC«IC*1 

DO  5i0  J»JqE.JPE 

M»»*l _ 

I F  <  H , LT . 4 )  GO  TO  501 

IC«IC*1 

M*1 


01  04-03-73  1C. 515 


501  jr< JCOUNTJJ) .EQ,0)  GO  TO  500 
Il»AND<  I  * ( I C  > • MASK ( M  > ) 
ll«ll/KONST(M) 

IF(  I1.GE.2048)  I1*0R(I1,KM) 

Jl«JCO(iNT(  J) 

Ll«LCOUNT{ J) 

I  DAT ( Ji .  1 1 )  *  1 1 
JCOJNT( J)«JCOUNT( J>»t 
500  CONTINUE 

000  continue 

NNP*4*NN 

NN4*NNP-3 

DO  900  KZsNN4,NNP,2 

I C«  IC*1 

KZP*KZ*1 

J0D(KZ  )  *  A  NO ( I A( IC) .MASkA )/IHALP 
900  JQD(KZP)«ANDUAUC)#MASKB) 

1000  I C* I C*5 

I T { I  PR  I  NT , EQ ■ 0 )  GO  TO  64 
URITE(6,45)(!0D(!).l«1.4) 

45  FORMAT(1HO,4(2X.010) ) 

kR I TE ( 6. 46 )( J0D< 1 >. 1*1,32) 

46  FORMATUOX. 012. 2110. 2X.  012) 

64  I JKs I JK*1 

CO  5001  1*1. NUM 

~  IFMjK.hE',NT!MES)  GO  jO  20 - 

J  JKiC 

REC0RD*RFC0RD*1, 

LO  40  Mtt.NUM 
FME AN*  0 , 

PO  6U  J*1 . 320 

- roATTrWTT7Rl - - - 

60  FMEAN*FDAT*FMEAN 
FMEAN=FMEAN/320, 

FAVG(M)*( (RECORD-1. )*FAVG(M)#FMEAN)/RECORD 
UR  I  TE ( 6 , 47 )  FHEAV.FAVG(M) 

47  F0RMAT(10X.14HRECGRD  MEAN  =  «Fi2,4.  1&H  PARTIAL  MEAN  « 

4g  iAvG(M)*f  AVG'(M)*n,5 - - 

20  IF<  IPARTY.nE.I)  GC  Tq  5001 
IND«KCOUNT(I ) 

URITE (6.311)  IND 
HI  FORMAT  (  5X ,  RHCHANMEL  .  I  ? ) 

KEND*JC0UNT( IND)-i 
TRITE (6. 40)  ( I  BAY  < J,  J) ,  J«l,KFND) 

5001  CONTINUE 
J  P ART  Y  »  0 
5000  RETURN 
ENC 

?T7^rwWT'cF  HtMftRV  USCD  BY  THIS  COMPILATION - - - - 
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«Br 


02  04-03-73  10.665 


1 

lbl 

SUBROUTINE  TO  RP  AD  BAHKI  TAPE 

2 

symdef 

BELLRD 

3 

SYHREF 

BSF 

oooooo 

4 

BLOCK 

JAZ  7 

oooooo 

5 

JZ 

BSS 

1859 

003503 

6 

IBORT 

BSS 

1 

0C3504 

7 

IPARTY 

BSS 

i 

OOOOOO 

8 

USE 

“PREVIOUS 

oooooo 

9 

BELLRD 

SAVE 

oooooo 

00(002710(00 

010 

000001 

00003563C030 

010 

000002 

000035754000 

010 

000003 

000035741000 

010 

000004 

777776  2350  07 

000 

10 

Tru 

■-S7DL 

000005 

000033  7550  00 

010 

11 

Sta 

StR 

000006 

000001  0010  00 

003 

12 

MME 

GE I  NOS 

000QQ7 

05  OCOO  uoooou 

000 

13 

rtb 

000010 

000031  000032 

311 

14 

ZERO 

FC » DCW 

0Q0Q33  QlOOOO 

010 

15 

ZERO 

STR 

000012 

0 0 0 0 L 2  OCIO  00 

30J 

16 

RME 

GEROAD 

000013 

000033  2350  00 

010 

17 

lda 

STR 

000014 

000040  3750  00 

010 

18 

ana 

,0070000000000 

000015 

OQ0041  1150  00 

010 

19 

CRPA 

<0030000000000 

000016 

000023  6000  00 

010 

20 

TZE 

adopt 

000017 

000042  1150  00 

010 

21 

CMPA 

*0040000000000 

000020 

000026  6000  00 

010 

22 

TZE 

EOF 

000021 

OOOOOO  2350  07 

000 

23 

LOA 

«0,OL 

000022 

000027  7100  00 

010 

24 

TRA 

RETURN 

END  OF  BINARY  CARD  SOBRQUT  J 

000023 

000001  2350  07 

000 

25 

ABORT 

lda 

*1 »  dl 

000024 

013504  7550  00 

030 

26 

STA 

IPARTY 

000025 

000030  7100  00 

010 

27 

TWA 

Rf7URN*1" 

000026 

OOOOCl  2350  07 

000 

28 

EOF 

LPA 

*1,0L 

000027 

013502  7550  00 

030 

29 

RETURN 

sta 

JZ*185B 

000030 

00000171(000 

010 

30 

RETURN 

BELLRD 

000031 

000000000001 

000 

31 

FC 

SC  I  • 

1,000001 

000032 

OlOOOO  0(3502 

030 

32 

DCW 

I  OTD 

J7.185B 

000033 

33 

STR 

BSS 

ERROR  LINKAGE 

000035  OOOOtOOOCOOO  000 

000036  222543435124 000 


LITERALS 


000040  070000000000  000 

080041  030000000000  OOU 

000042  O4000Q00CCOQ  OOP 

END  OF  BINARY  CARD  SUBROUT! 
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43  IS  THE  NEXT  AVAILABLE  LOCATION. 


END 
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PAGE 


1 

LflL  SUBROUTINE  TO 

WRITE  OUTPUT  TAPE 

(  7  TRACK  VERSION  ) 

2 

SYHDEF 

PROC, WEF, RSF 

oooooo 

3 

BLOCK 

JIM 

000000 

4 

JZ 

bss 

326 

oooooo 

5 

USE 

PREVIOUS 

oooooo 

6 

PROC 

SAVE 

0.1 

oooooo 

010 

000001 

000000220303 

000 

000002 

000000221003 

000 

000003 

00G333630G30 

010 

000Q04 

000333754000 

010 

000005 

000333741000 

010 

000006 

O00C01740U00 

010 

000007 

0  0  0  0  C  2741000 

013 

000010 

OOOOOO  2200  03 

000 

7 

ldxo 

•  0 ,  DU 

000011 

OOOOOO  2210  03 

000 

e 

BIG 

lPxi 

«0,DU 

000012 

010000  2360  10 

030 

9 

LOOP 

ldq 

JZ,0 

000013 

000015  6050  00 

010 

10 

TPl 

G01 

000014 

300001  1760  07 

000 

11 

SBQ 

•  1 ,  Dl 

0000l5 

000024  7360  00 

0  0  0 

12 

001 

"QLS 

2TT 

000016 

000022  7370  00 

000 

13 

LLS 

18 

000017 

000001  0600  03 

000 

14 

ADXQ 

si, DU 

000020 

010000  2360  10 

030 

15 

LOQ 

J7.0 

000021 

000023  6050  00 

010 

16 

TPL 

G02 

000022 

000001  1760  07 

000 

17 

SBQ 

si ,  DL 

END  OF  BINARY  CARD  SLBROUTI 

000023 

000024  7360  00 

000 

10 

G02 

OLS 

20 

000024 

000022  7370  00 

000 

19 

LLS 

18 

000025 

000070  7550  11 

010 

20 

sta 

I  A  ,  1 

000026 

000001  0600  03 

000 

21 

ADXO 

si, DU 

000027 

000001  0610  03 

000 

22 

ADX1 

«1.DU 

000030 

000243  1010  03 

000 

23 

CMPX1 

■1 63 , DU 

000031 

000012  6gl0  00 

010 

24 

TNZ 

LOOP 

000032 

000001  0010  00 

000 

25 

MME 

gejnos 

OOOQJ3 

15  0000  oooooo 

ooo 

26 

000034 

000064  000065 

Oil 

27 

FC , DOW 

000035 

000066  OOOOOO 

010 

28 

zero 

STR 

000036 

000002  0010  00 

000 

29 

MME 

geroad 

000037 

000001710000 

010 

30 

RETURN 

PROC 

000040 

31 

WEF 

SAVE 

000040 

000042710000 

010 

000041 

00033363C300 

010 

000042 

000333754000 

010 

000043 

000333741000 

010 

000044 

000001  0010  00 

000 

32 

MME 

geinos 

000045 

55  0000  020001 

000 

33 

WEE 

END  OF  BINARY  CARD  3UBR0UT I 

000046 

000064  OOOOOO 

Oio 

34 

ZERO 

7C 

000047 

000066  OOOOOO 

010 

35 

ZERO 

STR 

000050 

000002  0010  00 

000 

36 

MME 

GEROAD 

000051 

000041710000 

010 

37 

RETURN 

WEF 

000052  -  38  BSr  S*VH 
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000052 

00o0547ic.r00 

010 

000053 

JO  0333631/)  00 

'13,0 

003054 

J00333754U00 

010 

000055 

J0033374J  000 

Oil) 

000056 

JOClUil  OUltJ  00 

000 

39 

MMfc 

gejnos 

000057 

47  OlOO  u 2 0 0 0 1 

oou 

40 

BSF 

0000*0 

■J0OQ64  ouoooo 

010 

41 

ZERO 

re 

003061 

0 0 u J <> 6  Oa:OQn 

010 

4? 

ZERO 

STR 

0  0  3  Q  *2 

JOOui  2  Ot  in  O'J 

000 

47 

MME 

geroap 

0030*3 

000053711000 

01  J 

44 

RETURN 

BSF 

OOT064 

TOOOLOuOi>OiT2 

000 

45  FC 

dCl 

1,000002 

000065 

J  0  C-  o  7  0  0  u  0  ?  4  3 

01O 

46  new 

IOTO 

I  A, 163 

000066 

47  STP 

BSS 

2 

000070 

4R  {A 

BSS 

163 

ERROR  LINKAGE 


000333  OOL'OIOOOU .'JO  000 

0003.^4  47M4623Lc?Q  OHU 

END  Of  BINARY  CAR'D  SI  tiP OUT  I 

49  END 

336  IS  THF  NEXT  AVAILABLE  LOCATION, 

GMAP  VERS  I  ON/ ASSET  BL  Y  DATES  JMPA  1  101  71/1  02971  JMPB  1101 7i/io297l 

THERE  WFRF  NO  WARN  I NIG  FLAGS  IN  THE  ABOVE  ASSEMBLY 
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PAGE 


1 _ LPL  SUBROUTINE  TQ  WRITE  OUTPUT  TAP EJ^  9  TRACK  VERS IJN  ) 


2 

SYMDEr 

PROC,wE"r,Bsr 

000000 

3 

BLOCK 

JIM 

ooouoo 

4  Ji 

BSS 

326 

, 

000000 

5 

USE 

PREVIOUS 

OGiiuno 

6  PRCC 

SAVE 

000000 

not  01.2711000 

010 

oonom 

000646631000 

010 

000002 

000646754000 

010 

000003 

000646741000 

010 

000004 

000000  2200  03 

OOd 

7 

lpxo 

*c,pu 

000005 

000000  2710  03 

000 

B  BIG 

inxt 

*p»0u 

000006 

01( 000  2360  10 

030 

9  LOOP 

LOO 

07.5 

000007 

000011  61  50  (JO 

010 

10 

TPL 

G<U 

000010 

000001  1760  07 

000 

11 

SBQ 

*1  *  PL 

000011 

000024  7360  00 

POO 

1?  G01 

OLS 

20 

000012 

000072  7370  00 

Cl  0  0 

13 

LLS 

IB 

000013 

OOtOll  0<CU  03 

000 

14 

AOXO 

si.  DU 

00OQ14 

010000  2360  10 

n  3  j 

15 

LPO 

JZ,0 

000015 

00(01 7  6t 5U  00 

010 

16 

TPL 

GO  2 

000016 

nOLOll  176''  07 

000 

17 

SBO 

«1»0L 

000017 

000024  73 60  (JO 

0  0  0 

IB  002 

QLS 

20 

000020 

000022  7370  Od 

0  00 

19 

LLS 

18 

000021 

000103  7550  11 

010 

20 

STa 

I  A  1 1 

000022 

000001  OfuU  03 

0  0  li 

21 

ADXO 

*i,nu 

END  OP  BINARY  CARD  SohnOI'T! 

000023 

000001  0610  03 

00  0 

22 

APX1 

*1 .  Dll 

00^024 

000243  1(10  03 

0  0  0 

23 

C^PXl 

*163*  PL) 

000025 

0000(6  6(10  CO 

010 

24 

TNZ 

LOOP 

000026 

0000 11  0(10  00 

000 

25 

MME 

QF  I  NOS 

000027 

15  0(00  (COOOG 

000 

26 

WTB 

000030 

000152  00C154 

on 

27 

ZERO 

re, new ~ 

000031 

00015  7  0 1 r f 0 n 

010 

2  B 

ZFRO 

STR 

*  000032 

J00002  0(1"  P 7 

000 

29 

gfroap 

000033 

000155  2350  CJ 

01  J 

30 

LPA 

DCR9 

000034 

000156  7550  00 

010 

31 

STa 

DCW91 

Q0n0?5 

JOOut  0  22  0  03 

00  J 

3? 

lpxo 

*0 ,  DU 

000036 

000 00 U  2210  03 

00  1 

3  3 

lpxi 

*0  1  DU 

000037 

U1000U  2?e 1  10 

03  J 

34  L  OOP 9 

LPQ 

J7 1  0 

000040 

000024  7360  OC 

0  0) 

35 

01  s 

20 

000041 

000  02  0  73 7 C  00 

0  0  0 

36 

LLS 

16 

00C042 

010001  2360  1C 

030 

37 

LPQ 

J?4l,0 

000043 

U 0 0  0 2  4  7360  00 

0  0  0 

35 

QLS 

20 

000044 

00(020  7370  00 

POO 

39 

LLS 

16 

OOOQ45 

010002  2360  10 

030 

43 

LPQ 

J?4?,0 

END  OP  BINARY  CARD  SIHRQUTI 

000046 

000024  7360  00 

00  0 

4 1 

QLS 

20 

000047 

0000(4  7370  CO 

0  0  0 

42 

LLS 

4 

000050 

000426  75  5  )  11 

Old 

43 

STA 

JA,  J 

000051 

00(014  73/u  OJ 

00" 

44 

LLS 

12 

000052 

01.00(3  2360  10 

03" 

45 

LPQ 

J7*3,n 

0  0 o  o 53 

J 0 (  9 2  4  736  J  (Id 

00O 

46 

QLS 

20 

56 
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000054 

001020  7370  00 

000 

47 

LLS 

16 

000055 

010014  276f)  ifl 

o  3  o 

4fl 

lHO 

JZ*4,0 

00"056 

100024  7369  00 

000 

49 

QLS 

20 

000057 

900010  7 3/9  09 

ooo 

50 

LLS 

8 

000060 

0  00  4c  7  7559  11 

01C 

51 

sta 

J  A  ♦  1  »  1 

000061 

0 0 C 019  7370  00 

000 

52 

LLS 

6 

0  0  0  0  5  2 

ill  (J  0 1 5  2360  10 

030 

53 

LOQ 

J7»5,0 

000063 

U00024  7369  00 

000 

54 

QLS 

20 

000064 

•JOOOiO  7370  0  0 

non 

55 

LI  S 

16 

000065 

OK  0 1  6  2369  1  0 

030 

56 

LOQ 

JZ*6,0 

000066 

J0CQ24  7369  00 

0011 

57 

QLS 

20 

000067 

000014  7379  OU 

ooo 

5« 

LLS 

1? 

000070 

J00430  755)  11 

01  n 

59 

STA 

JA*2, 1 

END  00  BINARY  CARD  Si  8”0l  T! 

00O071 

900004  73/j  00 

ooo 

60 

LLS 

4 

000072 

0 1 0 0  u  7  2369  10 

030 

61 

lOQ 

JZ*7,0 

00O073 

0010*4  7369  09 

OOO 

62 

QLS 

20 

0  0  0  0  74 

9000*0  7379  00 

OOO 

6  3 

LLS 

16 

000075 

niOOU  2369  10 

03) 

64 

LOQ 

JZ»P,0 

000076 

JOl 0*4  7369  00 

0  0  (j 

65 

QLS 

20 

000077 

J0G02  9  7 si 7 no 

OOO 

66 

LLS 

16 

uoomo 

000431  755 9  11 

OKI 

67 

STA 

JA*3,1 

oonioi 

0 01  011  0600  03 

OOO 

6  fi 

A  0x0 

s9,0U 

0001P2 

J  0  0  5  0  4  H'C"  03 

OOO 

69 

C^PYO 

*324 , DU 

000103 

000106  6000  00 

010 

70 

T7fc 

OUTPUT 

000104 

0CC004  0610  03 

009 

7l 

A  D  X  1 

*4,pU 

000105 

0 C 0 0 3 7  71  09  PC 

010 

72 

T»A 

L00P9 

000106 

0 0 C 0 1  1  OllJ  fJ u 

ooo 

73  OUTPUT  M^E 

GF  I  NOS 

000107 

15  00  00  0  C  J 0 0 0 

OOP 

74 

WTq 

000110 

900153  000156 

on 

75 

ZpRO 

FC9.DCR91 

000111 

oooiti  ooonuo 

01P 

76 

ZERO 

STR9 

OUfl  112 

OOOUO 2  Oll'J  00  000 

77 

gfroatt 

000113 

occoci  ouo  oo 

ooo 

70 

M^E 

GF  I  NOS 

END  OF  BINARY  CAK'J  SIRROUTI 

000114 

55  0000  02PC01 

090 

79 

wFF 

00OH5 

000152  OCUCP 

010 

8 19 

ZFRO 

FC 

000116 

000157  01 0900 

OU 

81 

ZERO 

STR 

000117 

U000l2  0019  00 

000 

32 

NME 

geroao 

000120 

U  0  C  0  0  1 7 1 0  0  9  G 

010 

83 

RETURN 

PROC 

900121 

84  wEF 

SAVE 

000121 

900123710900 

Oil' 

000122 

0 01646631  COO 

010 

000123 

000646754000 

010 

000124 

000646741000 

mo 

000125 

OOCOll  0C1O  CO 

ooo 

65 

MMfc 

gF  I  NOS 

000126 

55  0000  (29001 

ooo 

86 

wFF 

000127 

QJG153  000901 

010 

87 

ZERO 

FC9 

000130 

OOCltl  QCOOOQ 

010 

BP 

ZFRO 

STR9 

000131 

'I0C0C2  00  10  00 

ooo 

89 

Mf'E 

GcROAD 

050132 

950122711090 

010 

90 

RF  TU R6J  Wf 

000133 

91  RSF 

S  A  V  F 

000133 

00013571091,0 

010 

57 
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00013.4 

aflL64663LMa 

OIL 

000135 

'l  0  C  6  4  6  7  5  4  ii  1 0 

010 

000136 

000646741000 

OIL1 

END  OF  BINARY  CARD  SLBRO'.'Tl 

000137 

OOGOLl  CG1Q  L 3 

000 

V? 

14  ME 

000140 

47  OLOU  (-20001 

OOP 

93 

B'E 

000141 

_.O.Q 412.2  JCWiiO _ 

010 

94 

ZERO 

000142 

000157  OOfU'OO 

11 0 

95 

ZERO 

000143 

U000G2  0(11  00 

ooc 

96 

mm 

000144 

OOCOuL  OLIO  UJ 

000 

97 

M”E 

000145 

47  0(30  121001 

000 

9* 

bSF 

000146 

J0( 153  OLC  '00 

010 

99 

ZFRO 

JIM14  7 

_Mkl61.  OLMOQ 

510 

ICO 

zero 

000150 

0 0 G 0 u 2  OLIO  (JO 

000 

1  Jl 

MMb 

000151 

00(1 3 471 L 300 

010 

10? 

UPTURN 

000152 

000000001002 

030 

10  3 

rc 

BC  l 

000153 

OOGOuOOOL  03 

00  0 

1(4 

rc? 

bCi 

000154 

10(163  010243 

010 

105 

oc^ 

i  cm 

Q0Q155.  MD.4i4.aL  022  (L 

QIC  , 

106 

QCW9 

I  o  r  d 

000156 

107 

ncw9i 

BSS 

000157 

10* 

STR 

dSS 

10H61 

109 

STR9 

bss 

000163 

no 

I  6 

BSS 

PC0476 

111 

JA 

yes 

error  linkage 

000646 

OOOOCIOOudO 

000 

000647 

475146232MQ 

000 

end  or  binary  card  su&routi 

_ _ _ _ 112 _  PI1  _ 

650  IS  THE  NEXT  AVAILABLE  LOCATION, 

GHAP  VERSION/ASSEI'BLY  DATES  JMHA  110171/102971 
THERE  WERF  NO  WA«f.  1  *'G  FLAGS  I N  T  HF  ABOVE  ASSEMBLY 


Gc I  NOS 

FC 

STR 

GFROAO 

SEJMS 

EC9 

STR9 

GERCAD 

b  of 

1,00000? 

1,000003 

1  A, 163 
J  A 1 1 4  A 
1 

2 
2 

163 
1 A  4 


JMPB  110171/102971 


ocioonnoooooo 


Q4-0U-73  17,796  SEISMIC  SENSOR  DATA  FEATURE  EXTRACTOR 

NINE  TRACK  MAGNETIC  TAfl  VERSION 

SEISMIC  SENSOR  DATA  FEATURE  EXTRACTOR 
NINE  TRACK  MAGNETIC  TAPE  VERSION 

ALBERT  H,  PROCTOR  12  SEPTEMBER  1972  HONEYWELL  635  FORTRAN  IV 
NCHAN  IS  THE  MAXIMUM  NUMBER  Or  DATA  CHANNELS  i'O  BE  PROCESSED 
PROM  ONE  RUN  OR  THE  NUMBER  OF  MAG  TAPE  PASSES, 

NRUNS  IS  THE  NUMBER  OF  RUNS  TO  BE  PROCESSED  IN  THE  NEXT  PASS, 

XTRESR  IS  THE  THRESHOLD  VALUE  roR  THE  ABSOLUTS  AVERAGE  DEVIATION 
FROM  the  mean  POR  EACH  1,024  SECOND  hxndom, 

IRON  IS  AN  ARRaI  CONTAINING  the  RUN  NUMBERS  TO  RE  PROCESSED  IN 
THE  NEXT  PASS, 

XCHAN  IS  AN  ARRAY  CONTAINING  THE  CHANNEL  NUMBERS  SO  BE  PROCESSED 
FOR  EACH  RUN  IN  THE  NEXT  PASS, 

DIMENSION  JHEAD(6)#XBUFt6144>» JlAIL ( 6 ) , ITEM P ( 0 ) , XCHAN  (  20  ) # X*UN ( 20 ) 
DIMENSION  X ( 6  1 44  )  | FEAT ( 100 ) 

EQUIVA1ENCE  (X,IBUF),  ( JHEAD(7),XBUF(  1) ) 

IFC"  1 

READ ( 5,  1001  )NCHAN 
D°  21  KCHAN*1,NCHAN 
READ(5,  1000)  NRUNS,  XTRESK 

1000  FORMAT ( 215 ) 

READ (3,  1001)  (IRUN(X),I«1, NRUNS) 

1001  FORMAT ( 4002 ) 

READ (5,  1003)  ( XCH AN C),Z>1, NRUNS) 

1003  FORMAT ( 4012 ) 

H RITE (6,1004)  NRUNS, ITRESK 
1004  FORMAT ( 1H0, 1X,7HRUNS  ■  , T3 , 3X , 8HITRESH  -,I3) 

WRITE? 6, 1008 )  { IRON ( I) ,1*1, NRUNS ) 

1006  FORMAT ( 1H0. 1X,4HRUNS,5X,  16(21,02) ) 
wRIXEfs, 1012)  (ICHAN(I),I«1, NRUNS) 

1012  F0RMAT(2X,8HCHANNELS,1X,16(2X,X2) ) 

CAIL  REu 
DO  20  1*1, NRUNS 

nvect-o 

1  CALL  RD9(JHEAD(  1 ) , IChAN  ( I ) , 0 ,0 » EFC , IE0F ) 

IF(IEOF,NE,0)  go  TO  1 
IF(JHEAD(3),EQ,IRUN(I))  GO  TO  2 
C*IL  FSF 
GO  TO  1 

2  ISTART.321 
JSTART-1 

3  ISHIFT-O 
DO  4 0  J-ISTART, 5120, 320 
lF(j,EQ,1)  GO  TO  31 
DO  30  11*1,6 
JII*J»II-7 

30  ITEMP(II)-IBJF{JII) 

31  CALL  RD9(IBUF(J-6).ICHAN(I),0,0,IFC,IEOF) 

IFJIEOF.NE.O)  GO  TO  20 
XFN,EQ,i)  GO  TO  40 

DO  U  11*1,6 
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JXI*J*II»7 
JTAXMII)»IBUF(  JII) 

U  lBUF{JXI)«ITEtlP(II) 

U  0  CORTXRUE 

DO  6  <?■  JST  ART  *  5120*  1024 

1-3*1023 

HBAH-JTAXKU) 

isuri-o 

DO  9  R-J,l 

5  ISUn»ISUM*XABs (XBUF(K)"HEAM) 
X3Un-XSUH/1024 

PRINT  2001, XSUK 
2001  FORnAT( 10X,20X5) 

IF{IS0M,GB,ITRESH)  GO  TO  6 
iSHiFT-L 

6  CORTIBUE 
JTAUtS)>JTAXM5)"4 
JTAUr6)«JTAU(6)-800 
IF(JTAIL(6) ,GE,0)  GO  TO  81 
jTAILf6)«1000+JTAIL(6) 
jTAlL(5)«JTAII,(5)"1 

81  pRXRT  1003,JTAXL,XCHAK(X) 

1005  FORHATJ 1X,012,2X,X5,2X,012,2X,4(I5,2X) ) 
X F f XSHXFT , EQ , 0 )  GO  TO  8 
JJ.1 

J3TART- 1 

XF(iSHIFT,GE,5120)  GO  TO  70 

lDXF»9120-ISHlFT 

jSTART-IDIF/1024 

J3TART«JSTART*102U+1 

lSHiFT-XDIF/320 

I3hXFT-5121-(ISHXFT*320) 

DO  7  7-l3HlFT,5120 

ibuf(JJ)»xbuf(J) 

7  JJ*0J*1 

70  l3TAtT"jj 

GO  TO  3 

8  DO  100  J-liSUO 
100  X(Jj»XBUF(J) 

1*1 

sum»o, 

SUH2-0, 

DO  10  J>1,  1024,512 

CAXt  RSQAR(X(J),X(b121),R) 

SUB1-*UH1*R 
10  SWH2-3UN2*R*R 
pBAT ( l ) «SUN 1/2 , 

1-1*1 

FEAT ( l ) -SUM2/2 • 

m*i 
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04-04.73  17,796  SIXSHXC  SEKSOI  DATA  FEATURE  BXJRACTOR 

NINE  HACK  NAOXETXC  TAM  VEXSXOl 

L»o  11  3*1025,5120,512 
C  All*  *SQAR(X(  J),X(5121),R) 

S«B1»XJ«1** 

11  SUH2*9UM2*R*R 
r*Ai^i)»j;UHl/io, 
l-t* 1 

F*AX(i;-SUM2/10. 

l-L*1 

1,2.^025 
DO  200  J-1,3 
H-?*(3/3) 

C*H  R0UBIE(X(L1),X(L2),H) 
li»tl*2048 
200  l2?L2*2Q6l 

DO  9  3-513,1024 

9  x(jjBtX(J-5l2)+X(J)+X(J+l536)*X(J*204i)+~{J+3664))/5( 

CAlt  HAR(1OV(X,6,U0i4,FEAT(l)  »FEAT(i*1  ) ) 

1-1*2 

CAtt  IZM>AT(X,2,rEA?U),PEAT(l«D) 

1-1*2 

CAll  TH*ESH(X,2,4,riAT(l) ) 
t»l*4 

CAll.  «PECVT(X,12,101,FEAT(L),FEAT(1*12)  ) 

1-1*13 

CA(l  PEKPXK(X,2,0,1,FEAT(l)) 
if*l*1 

CAIL  BAX«OP{X(513) , 6 , 40 , 4 , FE ATU ) , FE AT ( l*  1 )  ) 
l»l*2 

CAll  BXRRAT ( X ( 5 13) ,2,FEAT( L ) , FEAT( l*  1 ) ) 

1-1*2 

CAlt  THIX3H(X(513) ,2  »4,FEAT (l ) ) 

1-1*4 

CAli  SPlcVT(X{313),12,101,PIAl(l),F*AT(l*12)) 

1-1*12 

CAlt  PEXPXX(X(513),2,0»1,FEAT(1)) 

PREDT  2000, (FEAT(ZXD), 1X0*1, l) 

2000  FO*BAX{10X,10*10,3) 

FOICX  1005,JTAXL,XCHAX(X) 

M-1 

DO  12  3*1, l, 6 

U-<>*5 

XF(l1,0I.l)  ll-l 

pUXCH  2002, ( FEAT (ZXD),ZXD-J, 11 ),H 
2002  f0XBAT{6E12,4,X8) 

12  H*H*1 

NX ECT*X VECT* 1 
jSf ART" 1 

jstAit*i 

60  TO  3 
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20  HI*I  2003'NVICT 

2003  FOlnAT(20H  TOTAL  VECTOI3  ■  ,15) 

21  CALL  RE* 

STO» 

KID 

><><><><><><><>o<><><><><><><><><><><><><><><><><><>< 

ZEOr  DOES  MOT  APPEAR  IN  READ,  DATA,  COMMON  OR  LEFT  OF 

><><><><><><>o<><><><><><><><><><><><><><><><><><><>< 


><><><><><><><}<><><><><><><><><><><>o<><><><><><><>< 

R  DOES  NOT  APPEAR  ZN  READ,  DATA,  COMMON  OR  LEFT  OF 

><><><><><>OUO<><><><><><><)<><><><><><><><><><>(>< 


23842  WORDS  OF  riENORY  USED  BY  TKZS  COMPILATION 


EWUAL3  (.) 


EUOALS  (>) 


01  04*04-73 


17,111 


SUBROUTINE  DOUBLE  (  xf  X,  Nun  ) 

C  THIS  subroutine  enables  the  calculation  Of  THE 
i  C  PONE*  SPECTRUHS  OP  TWO  REAL  WAVEFORM  SEGflINTN 

C  NXTH  ONE  CALL  TO  THE  FPT  SUBROUTINE, 

DlnKNSXOM  X (  1024 >, Y {  1024 ), SI (  1024 ),S2(  1024 ),TAB( 780) 
XFfNUM-1)  30,25,30 
25  DO  27  0" 1 , 1024 

27  T<J)»0, 

30  CALL  FFT(  10,"1.0,X,r,Sl,S2,TAB) 

•  XF(NUM-I)  285,400,285 

285  DO  300  K»2 1 5 12 

l»1028-K 

i  A'XlKl+XU) 

B»X{K)-X<L) 

C»T  < KT*T <  Z. ) 

DMIHl'Ki) 

XlK)-?A*A«-D*D)*,25 
300  rlKi»(C*C*B*B)*,25 

xni-o. 

H5131-0. 

CO  320  K-2,512 
H*N*5 12 

i  320  X (H ) »T ( K ) 

RETURN 

.  400  DO  450  K»2,512 

A*X(K)*X(K) 

B*t { K ) *T ( K ) 
l  450  XlK|*A*B 

'  xn)-o, 

i  return 

I  END' 

XXXXXXXXXXXXXXXXXXXXXXXXXXX 

SI  DOES  NOT  APPEAR  IN  READ,  DATA,  CONNOR  OR  LEFT  OF  EUUA13  (■) 

:xxxxxxxxxxxxxxxxxxxxxxxxxxx 


iXXXXXXXXXXXXXXXXXXXXXXXXXXX 

S?  DOES  NOT  APPEAR  IN  READ,  DATA,  COMMON  OR  LEFT  OF  EUUALS  (») 

:xxxxxxxxxxxxxxxxxxxxxxxxxxx 


IXXXXXXXXXXXXXXXXXXXXXXXXXXX 

TAB  DOES  NOT  appear  in  read,  data,  COMMON  OR  LEFT  OF  EUUALS  (*) 

;xxxxxxxxxxxxxxxxxxxxxxxxxxx 
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SUBROUTINE  FFT  (  NST  A  G  E*  SION  ,  XR,  XI ,  SCR  AT  1 ,  SCR  AT2#  SCR  AT  3  ) 

C  ALBERT  H#  PROCTOR  21  JANUARY  1972  HONEYWELL  635  FORTRAN  IV 

C  NSTAGE  IS  THE  LOG  BASE  2  OF  N  WHERE  N  IS  THE  NUMBER  OF  DATA  POINTS 

C  TO  Be  PROCESSED, 

C  "SIGN  IS  TH'E'f RAWS FORM/IN VERSE  TRANSFORM  FLAG.’ 

C  SIGN  IS  -1,  FOR  THE  TRANSFORM  AND  1,  FOR  THE  INVERSE  TRANSFORM. 

C  XR  WILL  CONTAIN  THE  READ  PART  OF  EITHER  THE  INPUT  OR  OUTPUT  DATA, 

C  Xi  WILL  CONTAIN  THE  IMAGINARY  PART  OF  EITHER  THE  INPUT  OR  OUTPUT  DATA, 

C  SCR  AT  1  AND  SCRAT2  ARE  SCRATCH  ARRAYS  OF  LENGTH  N, 

C  SCRAT3  CONTAINS  THE  COSINE  TABLE  OF  LENGTH  3/4  N, 

DIMENSION  XR(2) ,Xi( 2) iSCRATI (2  5  »SCRAT2(2),SCRAT3(2) 

DATA  LSTAGE/O/ 

IF(SIGF)  12,11, 1 1 

11  ASSIGN  6  TO  ISIGH 
GO  TO  13 

12  ASSIGN  7  TO  ISIGN 

13  IF(NSTAGE-lSTAGE) 14,5*14 
14  lstage=nstage 

K"2**NSTAGF 
N2-N/2 
FLTNaN 

PHI2N»6.2831853/FLTN 
NPl»N2+1 
NPI 1 »NPI ♦ 1 
N4»N/4 
NPT2=N4*1 
N3FI2»3*N4-*1 
SCRAT3(  15*1. 

-  “ScRJfT3TrpT2T^C.  ' 

SCRAT3(NPI)«-1, 

SCR AT3 ( N3pl2 ) *0 , 

DO  1  I ■ 2 , N  4 
FI-I-1 

T5MP*ri*PHI2N 

tEmf-Cos  (tempi 

SCPAT3(I)*TEmP 
ISUB-N2+I 
ISUB1«NPH-I 
SCFAT3(ISUB)*-TEMP 
1  SCPAT3(ISUB1  )*-TEMP 

5  Lx  "1 

DO  3  J«1, NSTAGE 
NI»L 
L*2*L 
N2J*N/L 
NP*N  2 J 

-P5-T  TiTTFT - 

IN2J»(T-1)*N2J 
TK2K=IN2J+1 
IN2JI»IN2K+N4 
W1=SCRAT3(IN2K) 
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31  Q4-Q1-73  1,8.926 _  _ 

GO  TO  ISIGN* (6.7) 

6  W2»-SCRAT3(IM2JI) 

GO  TO  8 

7  W2»  SCRAT3(IM2JI) 

_ 8  DO  2  IR«1,NR _  _ 

isuB-nnnj 

ISUB1»ISUB+IM2J 
ISUb2*IS0b 1+N2J 
ISUB3«:SUfi+N2 

W8«M1*X8ilSUB2)-W2*Xl<ISUB2) 
wI"Lw2- XRJlsu?2J  1  * Xl ( ISUB2) 

SCR  AT  iTlSUB)  ■X~FTisUBlT*WR 
SCRAT2(ISUB)-XI(ISUB1)+WI 
SCRAT1(ISUB35-XR(ISUB1)-WR 
SCRAT2(ISUB3)«X:(ISUB1 )-WI 

2  CONTINUE 

DO  3  IR«1,N 
XR(IR)»SCRAT1(IR) 

3  XI(IR)«SCRAT2(IR) 
ir(SlGN) 10.9,9 

9  DO  4  IR«1,M 
XR(lR)»XR(lR)/FLTN 

_ 4_  XI(IR)«XI(  Ifll/FlTN _ 

10  RETURN 
end 

23647  WORDS  OF  MEMORY  USED  BY  THIS  COMPILATION 
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SUBROUTINE  HARMON  ( PS, NPSAKS , ISTART, MARGIN, $PACE#COUNT ) 
THIS  SUBROUTINE  DETERMINES  THE  MOST  FREQUENTLY  OCCURRING 
PAIRWISE  SPACING  BETWEEN  THE  LARGEST  PEAKS  IN  THE  POWER 
SPECTRUM  AND  THE  NUMBER  Of  TjMES  jT  OCCURRED, 

DIMENSION  PS  ( 5 12 ) , KCOUNT {  2 O')  ,  ITAB  {  20T#  IPREQt  iOT 

DO  9  I*  1  ( NPE AK3 

lpeak»istaRt 

PEAK*PS( ISTART1 

jstart«istart*i 

DO  1  J«JSTaRT»512 

IF(PEAK,3E»PS(J) ) GO  TO  1 

LPEaK»J 

PEAK«PS( j) 

continue 

IF  ( PEAK  ,  LT  ,  O'. )  GO  TO  17 

IFREQ(I)«LPEaK _ 

PS(LPIAK)«-P*AK 
JSTART-LPEAK+1 
DO  4  J*JSTABT,512 
PSJaPS(J) 

IF ( PS J , LT , 0  , )  30  TO  5 
0IF«ABS(PS(J-11 J-PSJ 
IF(DIF,GE,0. )  30  TO  4 
TEST=50,*(PEAK"ABS(P3<J-1) ) ) 

IF ( TEST , 3T « PEAK )  GO  TO  5 

PS(4)«-PSJ 

JJJ»j*1 

DO  3  JJ«  JJ  J,  512 _ _ 

PSJ»PS(JJ) 

IF(PSJ.LT.O. )  10  TO  5 
D I F  ■  A  B  S  (  P  S  ( J  J  ■  1 )  )  *  P  S  J 
IF(DIF)3^33,33 
PS  QJ  ji-PSJ 

GO  TO  5 _ 

PS<J )»-PSJ 
GO  TO  5 
JSfARI-JO+1 
PS( J)«-PSJ 
GO  TO  2 

JSTART«LPEAK-1 
DO  0  JJ«1,JSTART 
J-JSTART-JJ+1 
PSJ»PSf J) 

IF(PSJ,LT.0. )  GO  TO  9 
DlF«PSJ-*B5TPSTJ+in 
IF(DIF,LE,0.)  GO  TO  0 


IF(TSST.GT.PEaK)  GO  TO 

P5(J>-P5J 

JJJ»J-1 

DO 
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J-JJ.J-JJJJ+1 

P3.1*!>S(  J) 

IF ( PS J ,LT, 0, )  GO  10  9 
DIF»?SJ-ABS(PS( J+1)  ) 

IF ( DIF ) 77, 77, 7 
PS(J)«-PSJ 
GO  TO  9 
PS( J)»-PSJ 
GO  TO  9 
-JSTART«J-1 
PS{ J)—PSJ 

go- Tors 

CONTINUE 
GO  TO  18 
NPCAKS»I-1 
DO  10  I*  1 , 5  1 2 
PS ( I ) »ABS ( PS ( I )  ) 

ctct  s oitu^ttfiteoTn mKST - 

K.O 

DO  11  1*1,20 
KCOUNT ( I ) "0 
N 1»NP£AKS- 1 
DO  15  1*1, VI 

jstXrt«i*i 

DO  15  J«JSTaBT,NPEaKS 
IDir=IFREQ(J1-IFBEQ(I) 

MIN«IDIF -MARGIN 
MAX»IDIF ♦MARGIN 
IF(K.EO.O)  GO  TO  13 

Bo  12  KR-1.-R - - 

IF(ITAB(KK) ,LE. MAX. AND.ITAB(KK)  .GE, MIN}  GO  TO  1« 
CONTINUE 

IF ( K . GE , 20 )  GO  TO  15 

KaK+1 

ITAB(K)-IDIF 


14  KC0UnT(KK)"KC0UNT(KK)*1 

15  CONTINUE 
IND*  1 

MAXMUM«KCOUNT( 1 ) 

_ D016I-2.K  _ 

IND-I 

MAXmiM-KCOGNTdl 

16  CONTINUE 
COUNT-XCOUNrflNC  1 
SPACE»ITAB(IND> 
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SUBROBTIN*  SORTUP  (HAT,*) 

C  THIS  SUBROUTINE  ZS  CALLED  IT  HARBOR  TO  ORDIR  THE 

C  R  LARBIST  PEAKS  ACCORDXRS  TO  TnEZR  AMPLITUDES , 

ozhRrszow  zrayoo 

Nl»N-1 
BO  2  Z«1( VI 
ZHD»Z 

lITTH»HAr(X) 

JSfARt»Z*i 
DO  1  JmJSTART#  n 
XFflRAt(J) ,OE, LITTLE)  SO  TO  1 
X,ZTTX|R«XRAY  ( J ) 

XND»J 

1  CORTZVUI 
XTEnPaZRAT ( Z ) 

ZRAt(X)>LXTTlI 

2  XRAX (X»D)»ZtEBP 
RETURN 

END 

23709  WORDS  O'  HEHORY  USED  BY  THXS  COMPXLATXON 
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C 

c 

c 

c 

c 


1 


3 


4 


SUBROUTINE  ThRESh  "(PS,  j  START, LEVELS, Pj ST ) 
this  Subroutine  compares  the  amplitudes  or 

ELEMENTS  IN  THE  POWER  SPECTRUM  AGAINST  A  NUMBER 
or  EQUALLY  SPACED  THRESHOLDS  and  RETURNS  the 
NUMBER  or  ELEMENTS  WHJCH  P ALL  BETWEEN  THESE 

THRESHOLDS. _ _ _ 

DIMENSION  PS(S12>,PIST(LEVELS) 

PEAK*PS ( I  STAR  T  ) 

small*peak 

JSTART* 1ST  ART^l 
DO  1  J«USTART,5t2 

SMALL«AMIN1(SMALL,PS(  J> )  _  _ 

PEAK«AMAXl<PFAK,Ps< J)) 

DI V  IDE* < PEAK- SMALL )/rLOAT( LEVELS) 

CO  3  J*l. LEVELS 
P  1ST ( J ) *0, 

DO  4  I  * ISTART, 512 

1ND«(Pr(  I)-sMALL>/DIv1DE*1.  _  __ 

IP ( I ND . GT , LEVELS )  IND«LEVELS 
P I  ST ( IND  >  *p I  ST  < IND>*1. 

RETURN 

end 
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SUBROUTINE  R5QaRTX7ScRaT,R> - 

C  THIS  SUBROUTINE  CALCULATES  R »  WHERE  R  Is  DEFINED 

AS  THE  RATIO  OF  THE  MAXIMUM  DEVIATION  FROM  THE 
MEAN  TO  THE  AVERAGE  DEVIATION  FROM  THE  MEAN  JN  A 
GIVEN  1/2  SECOND  WINDOW  OF  THE  '  l ME  WAVEFORM, 

D I  MENS  1  ON  X( 512 ) »  SCRAT (512 ) 

SUM«0.  "  ' 

DO  10  l  *1 #  512 
10  SUM*SUM*X  < ! ) 

SUM*SUH/512, 

DO  20  1*1.512 
20  SCRAT  (  I  )«ARS(X(  I) -SUM) 

xMas«o .  - 

SUM*0, 

DO  40  L* 1 » 512 
XMAS* AMAX1( XMAS, SCR AT ( L) ) 

40  SUM*SUM*$CRAT(L) 

_ RjJ-512  .*XM*S)/SUM _ 

return 

END 
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SUBROUTINE  B FN'RaT  <T»S  .TsTaWT  ,TeaT1  ,FEA"T2 ) 

THIS  SUBROUTINE  CALCULATES  THE  RATIOS  OF  ENErSY 
BETWEEN  OME  AND  TWENTY  HERTZ  ANB  BETWEEN  FORTY-ONE 
AND  SIXTY  HERTZ  TO  THE  ENERGY  6ETWEEN  TWENTY-ONE 
AND  FORTY  HERTZ. 

DIMENSION  PS (512  ) _ 

SUMx*0, 

I  SUM 2» 0  , 

>  SUm3«0 . 

)  IST0P*ISTART*1R 

DO  1  I « I  ST  ART , 1  STOP 

!  _  __  SUM1,SUM1*PS< I )  _ 

I  SUM2«SUM2*PS(  1*20 ’) 

\  1  SUM3*SUM3*PS( 1*40) 

*  FEATlxSUMl/SUM2 

>  FEAT2sSUM3/SUM2 

’  return 

I _ _ _ _ 
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SUBNUUTlNt  a"tCVT\ra,IMVtC!  •  I  Ini  VCC  I  »KfcA  rtB) 

C  THIS  SUBROUTINE  CALCULATES  THE  RATIO  OF  ENERGY 

C  ABOVE  A  SPECIFIED  FREQUENCY  TO  THE  ENERGY  BELOW 

C  AND  ALSO  CREATES  A  VECTOR  WHOSE  COMPONENTS 

C  CONTAIN  NORMALIZED  ENERGY  VALUES  FRON  CONSECUTIVE 

C  FREQUENCY  BINS  OF  THE  POWER  SPECTRUM, 

DIMENSION  PS<512>»VECTINVEcT> 

DO  l  I * NVECT 
1  VECT ( I  )  *0 . 

SUM«0  . 

DO  5  1  ■! » NVECT 


J  *  5  *  1-  3 


2 

5 

6 

K«5*m 

DO  2  L* JiK 

VECT ( I >*VECT< I ) ♦PS ( L ) 

SUM*SUM*VECT( I ) 

DO  6  I  »1 » NVECT 

VEC T ( I  ) *VEcT ( I )/SUM 

- 

RRV3T«0, 

rgyat*o, 

DO  3  1*1#  ITH 

3 

RGYBTpRQYgT*PS( I) 
j*ITH^i 

DO  4  I.J.512 

.  4 

RGYAT«RGYAT*PS(  I  5 

REATEB»R6yAT/RGyBT 

RETURN 

ENP 

23715  WORDS 
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SUBROUTINE  PEKP1K(PS, IST4RT, THRESH, COUNT) 

THI*,  SUR«0UTlNfc  CALCULATES  THP  NyMBER  OP  PEAKs 
IN  A  POWER  SPECTRUM  WHICH  ARE  ABOVE  SOME  SPECIFIED 
PERCENTAGE  OF  THE  MAXIMUM  PEAK, 

DIMENSION  PS(5l2 ) 

XM 1  N*PS  (  I  START  ) _ _ _ 

AMAS«XMIN 
I N  J  T*  I  ST  ART *1 
DO  1  I  *  I N I T #512 
XMJNbAMINKXMIN.PSI  I ) ) 

1  XMASeAMAXl(XMAS,PS( I ) ) 

TEST«(XMAS-XM1N)«THRESH 
COUNT tO, 

PSMIN*PS( I ) — X  H I N 
IF(PSMIN.GT.TEST)  GO  TO  2 
JSW  ITsO 
GO  TO  3 

2  _ I  Sw  1 T « 1  _  _  _ 

3  DO  5  1*2,512 
PSMINsPSU  )-XMIN 

IF( ISWIT.EO.l)  GO  TO  4 
IF(PSmIN.GT.TEST)  ISWlTtl 
GO  TO  5 

4  1 F ( PSM I N . GT , TEST )  GO  To  5 _ _ _ 

I SW I T  *  0 

COUNT  *COUNT  *1 , 

5  CONTINUE 
RETURN 
END 
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1  LSI  TST  *IM,  SUBROUTINE  TO  SKI*  RECORDS  AND  REWIND  HAS  TAPES 


2  *  AL 

PROCTOR*  XSCP* 

RADC 

3 

5ia?z? 

rsr.Riw 

000000 

4 

rsr 

SAVE 

oooooo 

000002710000 

010 

000001 

000033630600 

010 

000002 

000033754000 

010 

000003 

OOOOJ3741VOO 

010 

000004 

000001  0010  00 

000 

5 

h;ie 

geznos 

000005 

45  0000  020001 

000 

6 

rsr 

000006 

000030  oooooo 

010 

7 

zero 

rc 

000007 

000031  oooooo 

010 

6 

zero 

STATUS 

000010 

000002  0010  00 

000 

9 

nm 

geroac 

000011 

000031  2350  00 

010 

10 

L00P1 

IDA 

STATUS 

000012 

000011  6050  00 

010 

11 

TPL 

LOOP  1 

000013 

000001710000 

010 

12 

RETURN 

rsr 

000014 

13 

REW 

SAVE 

000014 

000016710000 

010 

000015 

000033630000 

010 

000016 

000033754000 

010 

000017 

000033741000 

010 

000020 

000001  0Q10  00 

000 

14 

nits 

GEIWOS 

000021 

70  0000  020001 

000 

15 

REk 

000022 

J00030  oooooo 

010 

16 

zero 

re 

erd  or  bznart  carp  rsroooo2 

000023 

000031  oooooo 

010 

17 

ZERO 

STATUS 

000024 

000002  0010  00 

ooo 

11 

HME 

GEROAC 

000026 

000031  2330  00 

010 

19 

O 

o 

•Ef 

KJ 

IDA 

STATUS 

000026 

000025  6050  00 

010 

20 

TPL 

LOOP2 

000027 

000015710900 

010 

21 

RETURN 

REW 

000030 

0009000C0001 

000 

22 

rc 

BCZ 

1,000001 

000031 

23 

STATUS 

BSS 

2 

ERROR  UNXAGE 


000033  000000000000  000 

000034  264220202020  000 

E*D  Or  BINART  CARD  rsr00003 

24  END 

36  IS  the  NEXT  AVAILABLE  location, 

SNAP  VIRSION/ASSBmBLX  DATES  JMPA  110171/102971  JMPB  110171/102971  Jl 

THREE  VERB  NO  WARNING  flags  in  the  above  ASSEMBLE 
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1  LBX,  8EAD9TRK, SUBROUTINE  TO  RJAD  9  TRACK  DECODED  TAPES 


2  * 

JIM  ROACH,  IScP* 

RAoc 

3 

SYI1DEF 

RD9 

oooooo 

4 

USE 

PREVIOUS 

oooooo 

5  RD9 

SAVE 

0, 1,2, 3. 4,5 

OOOOJO 

000010710000 

010 

000001 

000000220003 

000 

000002 

000000221003 

000 

000003 

000000222003 

000 

000004 

000000223003 

ooo 

000005 

000000224003 

000 

00000* 

000000225003 

00) 

000007 

003045630000 

010 

000010 

003045754000 

010 

oooon 

003045741000 

010 

000012 

00000174)000 

01) 

000013 

000002741000 

010 

000014 

000003742000 

01  ) 

000015 

000004743000 

010 

000016 

000005744000 

010 

000017 

000006745000 

010 

000020 

000002  2350  11 

000 

6 

DA 

2.1 

000021 

003044  7550  00 

010 

7 

ST  A 

Bur 

000022 

000007  4500  31 

000 

8 

5TZ 

7,1* 

END  or  BINARY  CARD  READ9TRK 

000023 

000006  2350  31 

000 

9 

IDA 

6,1* 

000024 

003043  7550  00 

010 

10 

STA 

FC 

000025 

OOOOOO  2200  03 

ooo 

11 

LDXO 

=  0  ,  DU 

000026 

003041  2350  00 

01  1 

12 

IDA 

DUMMY 

000027 

000134  75b0  00 

01  ) 

T3 

STA 

TAUT 

000030 

003040  2350  00 

010 

14 

IDA 

DUMMY3 

000031 

000133  7550  00 

010 

15 

STA 

TALLYtJ 

000032 

000001  0)10  00 

00  ) 

16 

M'lE 

GEINOS 

000033 

0.1  0000  OOOOOO 

ooo 

17 

RTS 

000034 

003043  0)0135 

Oil 

18 

ZERO 

FC.DCW 

CRJW35 — OT3T35"0  57)000'  " 

TT1 

- 

'"ZERO  • 

STR 

000036 

000002  0)10  00 

00  J 

20 

MME 

GEROAn 

000037 

000136  2350  00 

010 

21 

IDA 

STR 

000040 

000002  7350  00 

00) 

22 

als 

2 

000041 

777400  3750  03 

ooo 

23 

ANA 

*0777400, DU 

000042 

211400  1150  03 

ooo 

24 

C'lPA 

•02 1 1400, )U 

- wo  oar 

-i 

°i 

o 

1 

TT) 

7  5" 

m  _ 

~IOOR 

000044 

000001  2350  07 

00) 

26 

IDA 

■1,DL 

000045 

000007  7550  31 

ooo 

27 

STA 

7.1* 

END  Of  BINARY  CARD  72AL9TRK 

000046 

000132  MOO  00 

01) 

23 

TRA 

EXIT 

000047 

003050  2360  00 

010 

29  100P 

LDQ 

■  1472 

WWW 

wo 006  7160  oo  ow 

W  ' 

QITS 

6 

000051 

000  134  7c.-20  06 

01  ) 

31 

STCQ 

TA^H ,  06 

000052 

003050  2360  00 

010 

32 

LDQ 

*1472 

000053 

000006  7360  00 

00) 

33 

gis 

6 

000054 

000133  7320  06 

01) 

30 

STCQ 

TAXlYg , 06 

"1 


<  EOF 


75 


2mi  02  Q4-Q3-73. _ 13,357 _ 


000055 

Ciy.  i  33 

2350  52 

010 

35  RUN 

IDA 

TALLYB, SC 

000056 

000134 

7550  52 

01  0 

36 

3TA 

TALLY, SC 

000057 

000055 

607)  00 

01  ) 

37 

TTF 

RUN 

000060 

000001 

0600  03 

00  ) 

38 

ADX0 

■  1 ,  DU 

000061 

000004 

1000  03 

00) 

39 

CMPX0 

■  4,  DU 

000062 

000047 

6010  00 

010 

40 

TNZ 

LOOP 

000063 

00)000 

2240  03 

00) 

41 

LDX4 

*0 ,  DU 

000064 

000000 

2230  03 

00  0 

42 

IDX3 

■  0,  DU 

000065 

000003 

236)  31 

00) 

43  START 

L00 

3,1* 

000066 

000124 

6000  00 

01 ) 

44 

TZE 

NEXT 

000067 

003051 

1760  00 

01  ) 

45 

S  B0 

■  1 

000070 

003052 

4020  00 

01  J 

46 

MPE 

■  109 

end  or  binary  card 

REAO9TRK 

000071 

000000 

6230  06 

00) 

47 

EAX3 

0.QL 

000072 

000000 

220C  03 

00  > 

48 

LDX0 

•  0,DU 

000073 

000140 

2360  13 

010 

49  BEGIN 

LD0 

I A  ,  3 

0  0  0  0  7  4 

003053 

2350  00 

01  ) 

50 

IDA 

■  0 

000075 

000022 

7370  00 

00) 

51 

US 

19 

000076 

003044 

7 5 5 C  70 

01  ) 

52 

STA 

buf.*o 

000077 

000001 

0600  03 

00  ) 

53 

ADX0 

■  1 ,  DU 

000100 

000022 

7370  00 

00) 

54 

LIS 

13 

000101 

000022 

7350  00 

00) 

55 

ALS 

18 

000102 

000030 

7310  00 

00) 

56 

ARS 

24 

000103 

003044 

7550  70 

01  ) 

57 

STA 

bjE, *3 

000104 

00000  1 

0630  03 

00) 

58 

ADX3 

■  1 ,  DU 

000105 

00000  1 

0600  03 

00) 

59 

ADX0 

■1,DU 

000106 

OOOOOO 

2220  03 

00  ) 

CO  LOOPX 

LDX2 

■  0 ,  DU 

000107 

003053 

2350  00 

01  ) 

61 

lda 

■  0 

000110 

00014C 

2360  13 

010 

62 

LOQ 

IA,3 

ooom 

000014 

T37'd  00 

000 

63  LOOPY 

"US 

~V2 

000112 

000030 

7350  00 

00) 

64 

ALS 

24 

000113 

000030 

7310  00 

000 

65 

ARS 

24 

end  or  binary  CARO 

READ9TRK 

000114 

003044 

7550  70 

010 

66 

STA 

buf,*o 

000115 

T00001 

0600  03 

00) 

67 

ADX0 

■  1  ,  DU 

000116 

0000  0 1 

062C  0l~ 

000 

68 

_ JTD3T2 

■  1  ,"DU 

000117 

000003 

1020  03 

003 

69 

CMPX2 

■  3, DU 

000120 

00C111 

6010  00 

01  J 

70 

TNZ 

LOOPY 

000121 

000001 

0630  03 

000 

71 

ADX3 

■  1  ,  DU 

000122 

000506 

1000  03 

00) 

72 

CMPX0 

■326, DU 

000123 

000106 

6010  00 

01  ) 

73 

TNZ 

LOOPX 

000124 

000001 

0610  03 

000 

Tu  NEXT 

ADX  1 

~  »T,T5U~ 

000125 

000506 

2250  03 

000 

75 

LDX5 

■126, DU 

000126 

003044 

0450  0 

01  ) 

76 

ASX5 

Bur 

000127 

000001 

0640  03 

000 

77 

ADX4 

*  1 ,  DU 

000130 

000003 

1 04f  03 

000 

78 

CMPX4 

■  3, DU 

000131 

000065 

60  ■*  C  00 

01  0 

79 

TNZ 

START 

000132 

00000171 0C00 

010" 

80  'tJTt 

“RETURN 

If&T 

coo  1 33 

81  TALLYB 

BSS 

1 

C  00 1 34 

82  TAUT 

BSS 

1 

000135 

o 

o 

o 

002700 

on 

83  DCW 

I0TD 

IA,  1472 

000136 

8  STR 

BSS 

2 

2Q915 

.02  04-03-73 

17.648 _ 

000140 

85  IA 

BSS 

1472 

003040  000140 

2700  40  010 

86  DUMKYfi 

tallyb 

IA,  1472,0 

END 

Or  BINARY  CARD 

READ9TRK 

003041  000140 

270C  00  010 

87  DUMMY 

TALLY 

IA, 1472,0 

003042 

88  LOCO 

BSS 

1 

C030T3- 

39  PC 

BSS ~  ' 

_  1 

003044 

90  BUF 

BSS 

1 

ERROR  LINKAGE 

003045  OOOOOOOOOOCO  000 

mmcr  siwiifinff  soo 

literals 

003050  000000002700  000 

003051  000000000001  000 

<103051  00'000<W3155 - TOO  ” 

003053  000000000000  000 

END  OF  BINARY  CARD  READ9TRK 

91  END 

3054  XS  THE  NEXT  AVAILABLE  LOCATION, 

GHAP  VERSION/ASSEMBLY  DATES  JMPA  110171/102971  JMPB  110171/102971 

THisi  whe  No  tf attstiT)"g  fi agt~itt t(Te  abovf  nnmr 
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o  o  o 


SKI  SHI C  SENSOR  FEATURE  TAPE  GENERATOR 

TUIS  IS  THi  CO:  lAoA  PRORRAM  WHICH  CONVERTS  DATA  CARDS  TO 

an  olra«s  oompatisle  magnetic  TAPE, 

PROGRAM  ALSOATA 

DIMENSION  X<l3o># IXIiSo* ,NUM(2o»«  INaME(j0),LFOR(2o>* IHEAD(7> 
EQUIVALENCES, IX) 

READ  1007*>SK1P 
1  o7  FORMAT! 1 9 1 

READ  1000*ITOR 
READ  looOiI^E 
ItOO  FORMAT ( 10 AB ) 

READ  1001SDJM,  JTOT,NUM 
NP3»NDIM*3  S  NLO»NDIM*i 
1  01  FORMAT  1 25 13) 

DO  7  I«1»IS<IP 
7  READ  TAPE  lg 

DO  l  IDUM"1* ITOT 
IT0P«NUM< IDJM) 

DO  2  I«l»  ITJP 
READ  loo««HEAD 

1„06  FORMAT(iX,U2,2X,I5,2X,012.2X.A<  IS.2X>> 

READ  LFOR*<X( J),J«i,N0IM) 

1^02  FORMAT<iS(lX,F3.0)) 

IX(A/3JM*l)«lHEADU)*100000*IMEAD«3)*UOOb*lHEAD(7>*1000*IHEAD(5) 

IX(NT'1M*2)»INAMEj  IDUM)  .AND.77B 

IX(N0IM*3)«INAME(IDUM> 

writ!  Tape  io# <xijj),jj«i,np3) 

2  CONTINUE 
1  CONTINUE 
PAUSE 

END  FILE  10 
REWIND  10 
PRINT  i003 
1  o’  FORMATIiHiI 

6  READ  TAPE  1o#(X(KK).kk>i,NP3> 

IF  (EOF  #io  M#5 

5  PRINT  1003# (X<KK),KK«l,NDRM) 

1  q3  FORMAT(iX,i0Ell,3) 

PRINT  100«*  (X(Kl'),KK-NL0,NP3) 

1  0R  FORMAT(lX#Ilo#2ClX,AI),//) 

GO  TO  6 
4  REWIND  io 
END 
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APPENDIX  B 


PROBABILITY  OF  CONFUSION  MEASURES 
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pair  c/h  rankjnjs 


RANKING 

1 

2 

3 

4 
3 
6 

7 

8 
0 

10 

11 

1? 

13 

1< 

13 

PAIR  C/A  RANKINGS 


measurement 

CONFUSION 

27 

.2l65970482 

3 

.3o6822612i 

,3732664996 

46 

29 

,388o8i3l44 

24 

,396491228i 

33 

•  40l0025()63 

7 

,4328599276 

43 

•  44()7i28933 

21 

•4484544695 

13 

,4613756614 

3 

,  49(j67j,i222 

4? 

,  5 o 3 2  o 245 06 

14 

•5l3li6i236 

44 

i5i4u86299 

39 

.5157894737 

ranking  MEASUREMENT  CONFUSION 


1 

47 

•I4g74n74g7 

2 

23 

,1777777778 

3 

46 

•1975589225 

4 

6 

•3o9175q842 

3 

33 

.3121212121 

6 

13 

.3153619529 

7 

24 

,3385942761 

8 

27 

,364983l650 

9 

28 

,38354377i0 

10 

42 

.3838383838 

11 

14 

,3832154882 

12 

36 

.3911616162 

13 

5 

.4061868687 

14 

39 

,407786i953 

13 

29 

.4292929293 
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pair  m/a  rankinss 


RANKING 

1 

2 

3 

4 

5 

6 
7 
I 
9 

10 

11 

12 

13 

1« 

19 

PAIR  M/T  RANKJN3S 


MEASUREMENT 

confusion 

4 

•0961497326 

3 

• o639q26738 

46 

•2376336898 

1 

•  29i67h23o 

45 

•265o4qio89 

2 

.309l977540 

24 

.3629679144 

23 

.3723262032 

31 

•3983957219 

3? 

•  4044u7647 

43 

•4i2433l95! 

44 

•4278074866 

22 

•47760695l8 

26 

•4799465241 

39 

•4986631016 

RANKfNG 

MEASUREMENT 

CONFUSION 

1 

43 

•02’46623o9 

2 

38 

• 0368i9i72i 

3 

44 

•0589324619 

4 

45 

• 08l3l8fl828 

5 

16 

•0958605664 

6 

42 

•096o7843l4 

7 

37 

•1326797386 

8 

22 

•1547930283 

•1552287582 

9 

21 

10 

29 

•183986928i 

11 

20 

•1*44226579 

12 

39 

•2139978213 

13 

23 

.2141612200 

1* 

15 

•22135q7625 

1’ 

7 

•243q827887 
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PAIR  T/A  RANKJN3S 


PAIR  M/H  RANK  1 N 3 S 


RANKING  MEASUREMENT 


1 

37 

2 

29 

3 

15 

4 

38 

5 

16 

6 

7 

7 

36 

8 

42 

9 

14 

10 

35 

11 

43 

12 

20 

15 

13 

14 

41 

15 

19 

RANKING 

MEASUREMENT 

1 

3 

2 

4 

3 

28 

4 

43 

5 

27 

6 

6 

7 

21 

8 

5 

9 

31 

10 

32 

11 

41 

12 

19 

13 

44 

14 

1 

13 

2 

CONFUSION 

.0’286i9529 

.1239898990 

•13484oo873 

.1654882155 

*2o”326599 

»21094276o9 

.245o3367qo 

*3o3o7239q6 

»3o7o286l95 

.3277356902 

.3739898990 

.39o74o74q7 

*4133838384 

.4287878788 

«4436026936 


CONFUSION 

•111455io84 

•1118209642 

.1258292791 

.I4i696i52i 

•1560703229 

.1620964175 

.l86808934i 

.2001879699 

.2449690402 

.24530o7519 

.26o67qo575 

•2757q76515 

.29o4i35338 

•3io482o876 

.34q5572755 
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PAIR  C/T  RANXJN3S 


PAIR  C/M  RAN* I N 3 S 


RANKfNG 

MEASUREMENT 

CONFUSION 

1 

38 

.03l85l85i9 

2 

18 

.0962962963 

3 

39 

.lo37fl37o37 

4 

17 

.2296296296 

3 

48 

.237o370370 

6 

37 

,2444444444 

7 

47 

,23i83i85iS 

8 

23 

,2666666667 

* 

49 

•274o74q74i 

10 

44 

.274074074i 

11 

24 

,3555535536 

12 

13 

,362962963q 

13 

23 

.3851831852 

14 

43 

•4000000000 

15 

2? 

.4222222222 

iNKfNG 

MEASUREMENT 

C INFUSION 

1 

47 

.14o2305447 

2 

43 

,2383877996 

3 

6 

.2655773420 

4 

29 

•2873638344 

5 

28 

.2956427013 

6 

26 

.3023965142 

7 

3i 

.3097494533 

8 

4 

.3177559913 

9 

32 

,3l79i939oo 

10 

3 

.3248910673 

11 

4  7 

.3321895425 

12 

26 

,346732q26i 

13 

27 

,37696o7843 

1* 

21 

,39o7932o7o 

13 

23 

. 398366ol3i 
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M 


PAIR  A/H  RANKIN3S 


RANK  fNQ 

measurement 

CONFUSION 

1 

5 

,1431989064 

2 

27 

»166q97o6o8 

3 

2* 

,2137303486 

4 

6 

,2252648667 

5 

«1 

,3320232399 

6 

19 

,34l379ol57 

7 

40 

,4427973342 

S 

1« 

,4636876282 

,46420o2734 

9 

43 

10 

21 

,4959415584 

,5037593985 

11 

47 

12 

17 

,5260167464 

13 

42 

,532o40328i 

i< 

3 

,5331937799 

,533706425i 

H' 

2  0 

pair  t/h  rankinss 


RANK  fNQ 

MEASUREMENT 

CONFUSION 

1 

16 

,1048175996 

.11?52i0248 

2 

38 

3 

29 

.1272626009 

4 

7 

,1490392648 

5 

37 

,1493734336 

6 

27 

,2309663046 

7 

15 

,26n8i87l34 

8 

5 

,2683375104 

9 

17 

,3206349200 

10 

39 

,3874129769 

11 

45 

»42oo5ol253 

12 

46 

,42534h3o6 

13 

1ft 

,44q9356725 

,4774i576i6 

14 

3 

19 

42 

,48o2flo5ol2 
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APPENDIX  C 

LINEAR  DISCRIMINANTS 
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***** 

FISHER  LOGIC 

***** 

*  **  ***  ** 

***** 

***** 

NODES  IN  SET 

Uel  h  arc  a  big  t  men  m 

l3l  C 

***** 

•AIR  1 

*  NODE  ARC  A  -  NODE 

HEL  h 

fisher  * 

coenciENTs 

3,968i3970E-o2  -9 . 68449i65E-04  -4 . S»3 i4 8 i E - 0 3  5,46598674p.Q2  -2,729217q3E-01 

6.973313296-03  -8 , 2«998958E-03  6.64553796E-02  2 ,90634536p-0i  -4 , 76369896E-02 

4,589982206-o2  6 . 04i66410E-0i  7.42977735E-q2  -4 , 7494 7i6iP-0i  -4,54l5269lE-oi 

l,5B537605E-oi 

THRESHOLDS 


3.9oo^^l35E  00 


***** 

RAIR  2 

♦  NODE  BIG  T 

-  NODE  HEL  H 

FISHER 

COEF I C I  ENTS 

4,697oB409E-o2  -1 . 93562543E-03  l ,4866g333E-o4  6 , 4462B8i8p-03  -4,76223T67E-02 

1.23o247i3E-01  -4 . 70777325E-03  3.53239344E-01  2.905i5473P-o1  5,562838i9E-o2 

1.q4*13262E-o2  -4 . t j)778l87E-oi  -5.2o48598oE-oi  “4.l0062826P-ol  -3.992o7l57E«*oi 
-2.792H913E-02 

THRESHOLDS 


•2.l7S37oooE  00 

•  PAIR  3 

*  node  big  t 

-  NODE  arc  A 

FISHER 

* 

* 

* 

***** 

2.847i3l67E-o2 
1 .i345qo48E-o1 
4,346o85i4E-o2 
4«6o8ol8846-o2 

4 .348673546-01 


6.0*5693986-03 
9 . 344i556qE-o  4 
3.7?8i5440E-0i 


COEFICIENTS 

1  ,o336i576E-o4 

6.3io709^6E-oi 

3.8oil5643E"oi 

THRESHOLDS 


3,B2463o38f-03 

4,86526098p-0l 

2.24040779p-0i 


8,76o89367e-o2 

1.89365433E-02 

7,54861956E-02 
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t 


FAIR  4 


*  node  men  m 


•  NODE  MEL  M 


F I  Smew 


COEF  JCIENTS 

2.3l4l278(jE-02  -8.  ?936u74E-o4  -4 .  B7289242E-o4 
-7,3558545*6-02  -2 . 3l964o98E-()3  -5.1142H64E-01 
■6«o3(j7543()E-02  4 , 75379822E-oi  2 . 853o5202E"oi 

"1 ,  o892l*44E*o2 

THRESHOLDS 

-9.642B5189E-01 


1.56999376F-Q2  2 1 456i,5n9ftE-oi 

-5.8363B325F.oi  -2,47784696e-o2 
5.5i895298f-02  1 , 28423779E-Q1 


♦  node  men  m 


node  ARC  k 


FISHER 


6,385g98i3E-o2 

•8»683Bioi56-o2 

-8,484o6n86E-o2 

4,7oBoB742E-o2 

2,357552676-01 


COEF I C I  ENTS 

■2 • 9?68889qE-o3  -lt22o36734E-o4 


1.63364582F-Q3 


■1 , $7395633E«o9  -3.6o53oo52E*oi  -3 , 926i24i6F-q1 
3«l54i3584E-oi  i.8o34B32lE“oi  4,559j)9276F-ol 

THRESHOLDS 


7.74779748E-02 

•4.5o*i373oE«o2 

5»890o7i28E-oi 


*  NODE  men  m  -  node  big  t 


fisher 


1.»203255’F-o3 

“1.789ii342F-o1 


COEF IC I  ENTS 

-5.5gB*4445E-02  5 , 7s9i n  0  ^  o  E • q  3  i.724fl67i 8E“e4  1.92o32559f-q3 

•1, 1B786829E-01  4,tfl996653E-fl4  -l.858l4o§5E-oi  -1.769n342F-ol 

9,757746B8E-o2  4 , 468{6622E-qi  4 ,8i252oi8E*oi  5,63346962P"o1 

1.57i87l28E-02 

THRESHOLDS 

1.4239o262E-oi 


3,15769992E-o2 

■2,15io6346E-o? 

3«9643oi46E-01 


♦  NODE  131  C 


-  node  hel  h 


FISHER 


COEFICTENTS 

9,8i44q973E-o2  -4 , 44773j.76E-03  i  4 688 1 n 5 5 E - q 4  -9 , 54345732F-()4  2»3B345"oi}E"01 
.  2«l8534455£«(^J  -9,Ai695737E-04  -2.62i43B02E-01  -7,8i897234f-o1  *1  ,lo®73436E-oi 
-2,935l6li46-01  2. 3iB8i2i6E-0t  1.00293l23E-0i  -1 , 48426786f-02  2,84237662E-0i 


9,54345732F-04  2.58345R00E-01 

7,8i897234f-o1  *1  »lo®73436E-oi 


1,525o218oE-02 

-1.45736931E-01 


THRESHOLDS 
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fai r  a 


♦  NODE  131  C 


NODE  *PC  A 


FISHER 


COEFICIENTS 

3.687ll7l0E-02  -1 , 93352649E-03  1 .11484727E-o4  -5 , 66834421*. e3 
2.6o489l4lE«o4  1 , 735-19o76E-o4  8 . 674895876-01  -4(34q33285f>o1 
-1.563879896-31  4,4B29i408E-02  4 , 92898694E-B2  7 , 954B2735f-02 


4.131191586-03 

1.778546196-01 


5,io8875loE-o2 
,5*3if>949o*E-o2 
1 »3B7B4B37E*oi 


THRESHOLDS 


FAIR  9 


♦  node  131  C 


•  node  big  t 


FISHER 


COEFICIENTS 

.29i24399E-o2  3 . 38280485E-o3  2 - 1 0 9 n 9*i 9E -0 4  -5,8i29967«F-B3 
i i6293q96E-o4  3.?i56?445E-b4  -6 . 88i6074iE-o1  -6 ,4B37i32oF-Ol 
,l4i54548E-oi  7 . 2B82o436E-B2  2.92250744E-02  6 ,  366!5539f-o2 
1 637oi439E-o3 


>637oi435E*o3 

i66589q73E-o2 


4i3ft942?B7E-B2 
2»87653i23E-Bi 
8  «6449j949E-B2 


THRESHOLDS 


FAIR  1B 


*  NODE  131  C 


NODE  MEN  H 


fisher 


COEFICIENTS 


-l*45i09285E-02  1 . l49i3285E-03  4 . H7l425lE-B5  -1 ,46447234F-02  3 . 46B32284E-o7 
2.9886336BE-o4  7 . Bi998367E-04  1 . 110124236-01  -7 . 586*2477p-0i  -2.9659ii45E-B1 
-5«602*83llE-01  -6.8R66i87lE-B2  -9 . o345o8676“B2  "3,53279527p-b2  4 . 21597498E-Q3 
2.277258186-03 

THRESHOLDS 


1,838595366-01 
END  of  THIS  LOGIC  SfT 
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